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FOREWORD 

This  final  report  describes  the  elements  and  events  occurring  in  a 
typical  ejection  with  a  state-of-the-art  escape  system,  the  theory  and 
formulation  of  the  simulation  model,  and  the  procedures  necessary  to 
successfully  prepare,  execute,  and  utilize  the  Generalized  Escape  System 
Simulation  (GESS)  Computer  Program  and  the  related  ACT  and  DRAS  programs.  The 
GESS  Programmer's  Manual,  a  companion  guide,  lists  the  annotated  FORTRAN-IV 
program  code  and  represents  the  second  of  two  volumes  of  GESS  documentation. 
A  portion  of  this  work  was  performed  by  Ketron,  Inc.,  in  accordance  with  NADC 
Contract  N62269-81-C-0206,  Task  No.  630-1944. 

This  work  was  sponsored  by  the  Naval  Air  Systems  Command  (AIR-5312)  under 


Airtask  No.  A5315312-001D-1531W1100,  Work  Unit  No.  WC530. 
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ABSTRACT 

The  Generalized  Escape  System  Simulation  (GESS)  program  is  a  comput¬ 
erized  mathematical  model  for  dynamically  simulating  the  performance  of  exist¬ 
ing  or  developmental  aircraft  ejection  seat  systems.  The  program  generates  six- 
dimensional  trajectory  predictions  of  the  aircraft,  seat/occupant,  occupant 
alone,  and  seat  alone  by  calculating  the  forces  and  moments  imposed  on  these 
elements  by  the  seat  catapults,  rails,  rockets,  stabilization,  and  recovery 
systems  included  in  most  escape  system  configurations.  User  options  are 
provided  to  simulate  the  performance  of  all  conventional  escape  system  designs 
undermost  environmental  conditions  and  aircraft  attitudes  or  trajectories. 
This  GESS  User’s  Guide  describes  the  elements  and  events  occurring  in  typical 
escape  systems,  the  theory  and  formulation  of  the  simulation  model,  and  the 
procedures  necessary  to  successfully  prepare,  execute,  and  utilize  this  and 
the  related  ACT  and  OKAS  programs.  The  GESS  Programmer's  Manual,  a  companion 
guide,  lists  the  annotated  FORTRAN  IV  program  code,  and  represents  the  second 
of  two  volumes  of  GESS  documentation. 
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1 .0  INTRODUCTION 

Designing  ejection  seat  systems  enabling  the  safe  escape  of  crew¬ 
members  from  high-speed  military  aircraft  is  a  formidable  task.  As  an  adjunct 
to  the  Navy's  on-going  escape  systems  program  of  system  engineering,  testing, 
and  incident  review,  a  mathematical  model  has  been  formulated  and  developed  to 
simulate  the  operation  of  any  ejection  seat-based  escape  system  under  most 
realistic  conditions  from  any  aircraft.  This  User's  Guide  desr  ,e.  the  pre¬ 
paration  and  use  of  the  Generalized  Escape  System  Simulation  .JS)  program. 

Originally  in  February  1972,  the  Naval  Aerospace  Recovery  Facility, 
El  Centro,  contracted  Vought  Aeronautics  Company  to  develop  a  simulation  model 
of  the  SIIIS-3  (now  SEU-3/A)  Escape  System,  made  by  Stencel  Aero  Engineering 
Corporation^^*.  This  program  was  a  reformulation  of  an  earlier  program  which 
Vought  produced  for  the  Air  Force  in  1970^).  This  Navy  program,  named  ICARUS^) 
provided  greater  detail  in  the  modeling  of  the  individual  system  elements  than 
did  its  predecessor.  These  enhancements  allowed  greater  flexibility  in  studying 
hardware  modifications  and  occupant  physiological  effects.  The  completed  pro¬ 
gram  provided  engineers,  civilian  management,  and  Navy  personnel  with  greater 
ability  to: 

.  establish  escape  system  performance  requirements, 

.  evaluate  competitive  escape  system  proposals, 

.  evaluate  escape  system  design  modifications, 

.  monitor  escape  system  test  programs, 

.  prepare  fleet  information  on  escape  systems,  and 
.  investigate  in-flight  escape  incidents. 


*  All  references  may  be  found  in  Section  7.0. 
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The  ICARUS  program  results  were  compared  with  test  data,  and  then  modified  to 
achieve  acceptable  results ^»5). 

In  1976,  the  ICARUS  program  was  obtained  by  the  Naval  Air  Development 
Center  (NADC) , Warminster ,  for  use  in  the  F-18  escape  system  development  program. 
This  effort  required  a  substantial  reprogramming  of  the  model  due  to  the  incor¬ 
poration  of  the  Martin-Baker  escape  system  in  the  F-18  aircraft.  As  this  diff¬ 
icult  task  continued,  it  became  clear  that  a  more  flexible  model,  which  could 
be  easily  configured  by  the  user  to  simulate  a  wide  variety  of  escape  systems 
and  aircraft  designs,  was  critically  needed. 

In  1980,  NADC  contracted  the  Computer  Sciences  Corporation  to  reformu¬ 
late  and  code  the  simulation  model  which  was  to  become  known  as  GESS.  This  pro¬ 
gram,  later  revised  and  documented  by  KETRON,  INC.,  provides  the  user  with  a 
greatly  enhanced  capability  to  describe  escape  system  parameters  without  code 
modification,  and  to  execute  simulation  runs  with  significantly  improved  compu¬ 
tational  efficiency.  The  capability  to  produce  graphics  output  through  the  use 
of  the  DRAS  program  has  also  been  added. 

This  GESS  User's  Guide  is  Volume  I  of  two  volumes  of  GESS  documenta¬ 
tion.  Volume  II  is  the  GESS  Programmer' s  Manual ,  which  includes  the  annotated 
GESS  program  listing. 
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2.0  ESCAPE  SYSTEM  OPERATION 

Following  an  aircraft  failure  requiring  emergency  egress,  a  typical 
sequence  of  escape  system  events  could  be: 

(1)  trigger  mechanism  initiation 

(2)  catapult(s)  ignition 

(3)  seat  first  motion 

(4)  drogue  gun/container  deployment 

(5)  catapult(s)  separation 

(6)  rail(s)  separation 

(7)  rocket(s)  ignition 

(8)  seat  separation  from  aircraft 

(9)  stabilization  system  operation 

(10)  rocket (s)  burnout 

(11)  recovery  chute  deployment 

(12)  recovery  chute  line  stretch 

(13)  seat/occupant  separation 

(14)  recovery  chute  full  inflation 

Although  every  supplier  of  escape  systems  incorporates  somewhat  different  con¬ 
figurations  and  mechanisms  for  accomplishing  these  tasks,  the  general  function¬ 
al  sequence  includes  most  or  all  of  these  events. 

The  components  of  a  typical  ejection  seat  system  are  shown  in  Figure 

2-1. 
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Upon  catapult  Ignition,  pressure  Is  rapidly  created  within  the  cata¬ 
pult  tubes  providing  thrust  to  Initiate  seat  motion  along  a  prescribed  path 
out  of  the  aircraft.  When  the  seat  headbox  enters  the  airstream,  after  either 
canopy  jettison  or  penetration,  the  drogue  parachute  is  deployed  from  the  seat. 
Prior  to  or  at  the  time  of  seat  separation  from  the  aircraft,  the  rochets  loca¬ 
ted  on  the  seat  are  ignited,  providing  additional  thrust  to  move  the  seat  and 
its  occupant  away  from  the  aircraft. 

As  the  seat  continues  its  trajectory  from  the  aircraft,  it  must  be 
stabilized  through  the  use  of  a  system  such  as  the  Directional  Automatic  Realign¬ 
ment  of  Trajectory  (DART)  system.  This  system  consists  of  a  brake  attached  to 
the  seat  through  which  one  or  two  aircraft-attached  lines  are  played  out,  pro¬ 
ducing  a  moment  designed  to  correct  any  adverse  seat  attitude.  These  lines 
break  when  the  seat  reaches  a  specific  distance  from  the  aircraft. 

Deployment  of  the  recovery  system  is  usually  accomplished  either  by  a 
drogue  parachute,  a  Wind  Oriented  Rocket  Deployment  (WORD)  system,  mortar  de¬ 
ployment,  or  a  combination  of  these.  The  time  of  deployment  of  the  recovery 
system  is  controlled  by  the  seat  sequencing  mechanism,  which  senses  seat  air¬ 
speed  and  altitude. 

As  the  parachute  deploys,  the  risers  and  suspension  lines  develop 
tension  from  the  differential  momentum  between  the  inflating  parachute  canopy 
and  the  seat /occupant.  This  causes  the  seat  and  occupant  to  separate.  The 
seat  moves  away  from  the  occupant  as  the  occupant/parachute  system  decelerates 
further  under  the  action  of  the  inflating  parachute  canopy. 
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3.0  FORMULATION 


3.1  Coordinate  Systems 

Trajectories  generated  by  the  Generalized  Escape  System  Simulation 
(GESS)  program  are  calculated  with  respect  to  specific,  "right  hand  rule," 
coordinate  systems  coinciding  with  center  of  gravities  or  other  specified  ref¬ 
erence  points  of  the  various  system  elements.  Each  of  the  program  coordinate 
systems  is  defined  by  six  (6)  degrees  of  freedom: 

A.  Linear : 

1.  X  -  forward  positive  displacement 

2.  Y  -  left  positive  displacement 

3.  Z  -  upward  positive  displacement 

B.  Angular : 

A.  Yaw  (R)  -  leftward  positive  rotation 

5.  Pitch  (Q)-  nose  downward  positive  rotation 

6.  Roll  (P)  -  right  wing  downward  positive  rotation 

These  coordinate  systems,  illustrated  in  Figure  3-1,  are  defined  as  follows: 

•  Earth-Fixed  Coordinate  System  (EFCS) 

The  EFCS  is  a  3-axis  orthogonal  coordinate  system  with  origin  at 
a  fixed  point  located  on  or  near  the  surface  of  the  earth.  Since  all  simu¬ 
lated  distances  are  relatively  small  compared  to  the  earth's  radius,  the 
errors  associated  with  neglecting  the  curvature  of  the  earth's  surface  are 
considered  negligible.  The  movements  and  rotations  of  all  simulated  system 
elements  can  be  described  with  respect  to  (wrt)  this  fixed  coordinate  system. 
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•  Aircraft  Coordinate  System  (ACS) 

The  ACS  Is  a  3-axis  orthogonal  coordinate  system  with  origin  at  a 
fixed  point,  usually  the  aircraft's  center  of  gravity.  All  linear  and  angular 
movements  of  the  aircraft  are  determined  by  establishing  the  position  of  this 
coordinate  system  wrt  the  EFCS.  The  position  of  all  other  system  elements, 
prior  to  and,  optionally,  after  the  separation  of  the  seat  from  the  aircraft, 
is  determined  wrt  to  the  ACS. 

•  Rail  Coordinate  System  (RCS) 

The  RCS  is  a  3-axls  orthogonal  coordinate  system  which  is  defined 
and  fixed  wrt  the  ACS.  Its  origin  coincides  with  the  mid-point  between  the 
rails  of  the  lateral  (Y)  axis  defined  by  the  two  lowest  rail  points  in  the 
ACS,  and  its  vertical  (Z)  axis  is  parallel  to  the  rails.  The  rail  system 
serves  to  restrict  the  motion  of  the  seat  to  a  specific  path  as  it  moves  out 
of  the  aircraft  cockpit. 

•  Seat  Coordinate  System  (SCS) 

The  SCS  is  a  3-axis  orthogonal  coordinate  system  with  origin  usually 
corresponding  to  the  initial  location  of  the  RCS  origin  wrt  the  seat.  The 
center  of  gravity  (CG)  of  either  the  seat/occupant  (S/0)  combination  or  the 
seat/alone  (S/A)  is  referenced  wrt  the  SCS,  allowing  the  trajectories  of  these 
system  elements  to  be  tracked  wrt  either  the  ACS  or  the  EFCS.  Points  of  appli¬ 
cation  of  the  various  forces  acting  on  the  S/0  or  S/A  are  also  referenced  wrt 
the  SCS,  allowing  the  calculation  of  moments  and  rotation  around  the  respective 
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•  Occupant  Alone  Coordinate  System  (OACS) 

The  OACS  is  a  3-axis  orthogonal  coordinate  system  with  origin  at 
the  center  of  gravity  of  the  seat  occupant  after  the  occupant  separates  from 
the  seat.  The  OACS  is  currently  defined  only  by  the  three  linear  degrees  of 
freedom  because  the  disjointed  non-rigidity  of  the  occupant  alone  makes  any 
angular  determination  difficult,  if  not  irrelevant.  Provisions  have  been  made 
for  the  future  incorporation  of  OACS  angular  degrees  of  freedom  if  such 
Information  is  necessary  at  a  later  date. 

•  Thrust  Vector  Control  Coordinate  System  (TVCCS) 

The  TVCCS  is  a  3-axis  orthogonal  coordinate  system  with  origin  at 
the  common  intersection  of  the  rocket  line  of  thrust  and  the  centerlines  of  its 
gimbals.  This  coordinate  system  establishes  the  orientation  of  the  rocket 
line  of  thrust,  which  is  modified  dynamically  as  part  of  the  vertical-seeking 
maneuver . 


Coordinate  systems  and  vector  directions  should  be  carefully  con¬ 
sidered  when  preparing  simulation  inputs.  All  output  reports  are  referenced 
accordingly. 
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3.1.2  Notation  and  Variable  Identification.  For  brevity,  vector  notation 
is  used  in  the  formulation  wherever  possible.  Individual  component  equations 
are  occasionally  used  for  clarity.  Position  vectors  are  denoted  by  the  symbol 
?.  Hence,  ?ASO  represents  the  position  vector  from  the  aircraft  CG  to  the 
seat/occupant  CG,  while  represents  the  position  vector  of  the  seat  alone 

with  respect  to  the  EFCS.  Direction  cosine  matrices  are  denoted  by  the  symbol 
j .  Hence,  D^g  represents  the  direction  cosine  matrix  which  transforms  compo¬ 
nents  of  a  vector  presented  in  the  AFCS  to  components  presented  in  the  EFCS. 

Several  naming  conventions  are  used  to  identify  variables  used  in  the 
formulation.  Refer  to  Section  4.0  for  specific  variable  descriptions.  These 
naming  conventions  are  summarized  below. 

SA  -  seat  alone 

OA  -  occupant  alone 

SO  -  seat/occupant 

NPTS  -  number  of  points 

RK  -  rocket 

WGHT  -  weight 

X  -  X  axis 

Y  -  Y  axis 

Z  -  Z  axis 

P  -  roll  rotation 

Q  -  pitch  rotation 

R  -  yaw  rotation 

VEL  -  velocity 

PORO  -  porosity 

recovery  parachute 


REC 
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DRO/DR 

-  drogue  parachute 

DRT 

-  DART 

POS 

-  position 

IGN 

-  ignition 

CAT 

-  catapult 

TVC 

thrust  vector  control 

CG 

-  center  of  gravity 

The  variable  names  used  in  che  formulation  have  been  selected  to  fa¬ 
cilitate  the  user's  understanding.  There  is  limited  correlation  between  these 
variable  names  and  the  coded  variable  names  used  in  the  actual  program. 

3.1.3.  Units  of  Measure.  Either  English  or  metric  units  of  measure  may  be 
used  in  the  GESS  program.  The  specific  units  required  for  either  system  of  mea¬ 
sure  selected  by  the  user  are  detailed  for  each  input  variable  in  Section  4.0. 
No  units  of  measure  are  referenced  in  the  formulation. 

3. 2  Vector  Mathematics 

3.2.1  Position  Vectors.  A  position  vector  is  a  representation  of  the  loca¬ 
tion  of  a  specified  external  point,  with  reference  to  a  stated  initial  point, 
characterized  by  both  distance  or  magnitude  and  direction  components.  By  stat¬ 
ing  the  location  of  a  point,  denoted  (2),  in  terms  of  the  origin  of  a  specified 
coordinate  axis  (l),  the  position  vector  (?12^  fro®  the  origin  of  the  coordi¬ 
nate  system  to  the  point  is  defined.  The  calculation  of  che  position  vector 
*23*  between  the  two  points,  (2)  and  (3),  is  facilitated  by  vector  subtrac¬ 
tion.  These  principles  are  illustrated  in  Figure  3-2. 
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The  definitions  of  the  position  vectors  used  in  GESS  are  given  in  con¬ 
junction  with  discussions  of  the  formulation  in  which  they  are  used. 

3.2.2  Direction  Cosine  Matrices.  Direction  cosine  matrices  (DCMs)  are  used 
to  transpose  position  vectors  from  one  coordinate  system  to  the  orientation  of 
another  coordinate  system.  Such  transformations  are  necessary  prior  to  vector 
addition  or  subtraction  operations  intended  to  determine  resultant  vectors. 
Hence,  the  addition  of  the  position  vectors  ?i2  ao<*  ^23  will  yield  the  new  posi¬ 
tion  vector  ?i3  after  ?23  has  been  transposed  into  the  ?i2  coordinate  system. 
Further,  given  position  vectors  fj2  end  £34,  the  vector  ?24  can  he  found  by  sub¬ 
tracting  ?i2  fro®  the  addition  of  Pj3  and  the  transposition  of  £34  into  the 
2  coordinate  system.  These  principles  are  illustrated  in  Figures  3-3  and  3-4. 
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The  generalized  DCM,  D21,  which  transforms  a  vector  in  the  initial 
coordinate  system  to  the  coordinates  of  an  external  coordinate  system  is  given 
by: 


dll 

d12 

d  13 

°21  ” 

d  2 1 

d22 

d23 

d31 

d32 

d33 

where: 


dn  *  cos  8  cos  y  (l.a) 

d  12  ■  cos  8  sin  y  ( 1 .b) 

d 13  “  -sin  8  ( 1  - c ) 

dzi  *  -sin  y  cos  a  +•  sin  a  sin  8  cos  y  (l.d) 

^22  *  cos  y  cos  a  +  s^n  a  sin  8  sin  y  (l.e) 

d23  a  cos  8  sin  a  (  1  .f ) 

d 3 1  =*  sin  a  sin  y  +  sin  8  cos  a  cos  Y  (l.g) 

d32  ■  -sin  a  cos  y  +  sin  8  cos  a  sin  y  (l.h) 

833  *  cos  a  cos  8  (  1 .  i ) 


and  where  the  modified  Euler  angles  are: 

a  =  positive  rotation  to  right  about  the  X"  axis  (Roll) 

3  =  positive  rotation  nose  downward  about  the  Y'  axis  (Pitch) 

Y  =  positive  rotation  to  left  about  the  Zi  axis  (Yaw) 

These  angles  are  illustrated  in  Figure  3-5. 
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Thus,  Che  position  vector  resulting  from  the  addition  of  the  position 
vector  to  the  external  coordinate  system  origin  to  a  position  vector  in  that 
external  coordinate  system,  as  in  Figure  3-3,  is  given  by: 

*13  *  *12  +  °21  *23 

And,  the  position  vector  found  by  subtracting  two  separately  coordinated  posi¬ 
tion  vectors,  as  in  Figure  3-4,  is  given  by: 

*24  ■  *13  +  °21  *34  -  *  12 

A  final  principle  involves  the  use  of  the  transpose  of  the  DCM,  which 
is  defined  as: 


dll 

d2i 

d31 

T 

d21“ 

d12 

d22 

d32  *  °12 

(2) 

d13 

d23 

d33 

If  the  DCM  which  transforms  the  first  coordinate  system  into  the  second  is  used, 
then  Figure  3-3  can  be  represented  as: 

*13  "  *23  +  °21  *12 

and  Figure  3-4  is  represented  as: 

*24  *  *31  +  °21  *12  “  *34 

This  is  possible  since  all  DCMs  are  orthogonal  matrices;  for  this  class  of 
T  —1  —1 

matrices,  *  D2J,  where:  D2J  =  inverse  of  D2^«  Thus,  a  single  DCM  is  used  to 
enable  position  vector  transformation  in  either  direction  between  two  separate 
coordinate  systems. 

The  definitions  of  the  specific  modified  Euler  angles  describing  the 
relative  orientations  of  the  various  program  coordinate  systems,  and  their 
respective  DCMs,  are  described  in  conjunction  with  the  discussion  of  the  formu¬ 
lation  in  which  they  are  used. 
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3.2.3  Quaternions .  Calculating  the  components  of  the  DCMs  for  each  of  the 
coordinate  system  pairs  using  the  modified  Euler  angles  involves  a  considerable 
use  of  trigonometric  functions.  The  cost  of  using  these  functions,  in  terms  of 
computer  time  and  CRJ  memory,  is  excessive  in  a  large  iterative  program.  This 
problem  has  been  substantially  resolved  in  GESS  through  the  use  of  unit 
quaternions. 

Quaternions  are  represented  symbolically  as  the  sum  of  a  vector  and  a 
scalar,  and  may  be  interpreted  geometrically  by  considering  them  as  operators  on 
vectors.  In  GESS,  unit  quaternions  are  used  to  update  DCMs  during  integration  of 
each  system  element's  equations  of  motion.  During  program  initiation,  quater¬ 
nion  components  are  initiated  for  each  DCM  as: 


*0  -  (1/2)I|<1  +  du  +  d22  +  d33)|]  1/2  (3. a) 

X1  *  <d23  -  d32>/^o  (3-b) 

>2  *  (d31  “  d13^/^o  (3.c) 

X3  *  (d12  "  d2l)/^o  (3.d) 


where  the  DCM  components  indicated  have  been  initialized  by  the  trigonometric 
calculations  used  in  equation  ( 1) . 

Following  initiation,  the  quaternion  is  updated  during  every  integra¬ 
tion  time  step  by  the  matrix  differential  equation: 
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where : 

p  s  angular  velocity  of  roll  angle,  a 
q  =  angular  velocity  of  pitch  angle,  8 
r  =  angular  velocity  of  yaw  angle,  y 


The  new  orientation  of  the  system  element  with  respect  to  the  DCM  reference 
coordinate  system  is  then  obtained  by  integrating: 

Xi  -  p^d t  (i  -  0, 1,2,3)  (5) 

And,  the  new  DCM,  D21  is  found  from: 


dll  <*12  d13 
°21  *  d2l  d22  d23 
d3l  d32  d33 


(6) 


where 


d'n  -  2(A^  +  aJ)-1  (6. a) 

dl2  *  2CX fX 2  +  1()*3)  (6.b) 

d13  “  2(1^13  -  AqA2)  (6.c) 

djij  m  2(A ^2  —  AQA3)  (6.d) 

d22  *  2<x0  +  X2)_1  (6.e) 

d“22  m  2(12X3  +  AgAi)  (6.f) 

d31  "  2(1^13  +  AqA2)  (6*g) 

d'32  -  2(A2A3  -  AqA  j)  (6.h) 

d33  -  2(Aq  +  x|)-l  (6.i) 
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Finally,  the  new  modified  Euler  angles  can  be  found  from: 


a'  ■  Can  1  (dJj/djj)  (7. a) 

S'  -  sin-1  (d|3)  (7.b) 

r  -  tan"1  ( 7.c) 

A  more  detailed  discussion  of  the  theory  and  use  of  quaternions  can 
be  found  In  most  linear  algebra  texts. 


Equations  of  Motion 


3.3.1  Aircraft,  the  aircraft  equations  of  motion  are  based  upon  linear  and 
angular  input  values  of  position,  velocity,  and  acceleration  for  each  of  six 
degrees  of  freedom  with  respect  to  Che  EFCS.  No  external  forces,  such  as  grav¬ 
ity,  thrust,  or  moments  are  considered,  although  the  linear  and  angular  accel¬ 
eration  input  can  be  structured  as  a  function  of  time  to  represent  nearly  any 
aircraft  maneuver.  The  input  values  are  then  integrated  over  the  simulation 
time  increment  to  determine  Che  linear  and  angular  time  history  of  the  aircraft 
with  respect  to  the  earth. 

The  aircraft  equations  of  motions  are: 


®AC  ■  /  sACdt 

(8. a) 

VAC  “  I  ^ACdt 

(8.b) 

®AC  "  /  ®ACdt 

( 8.c) 

“AC  “  /  “ACdt 

( 8.d) 
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linear  displacement  vector  wrt  EFCS 
linear  velocity  vector  wrt  EFCS 
angular  displacement  vector  wrt  EFCS 
angular  velocity  vector  wrt  EFCS 

3.3.2  Seat/Occupant  and  Seat  Alone.  Both  the  seat/occupant  and  seat  alone 
configurations  are  considered  as  rigid  bodies.  Consequently,  the  equations  of 
motion  for  both  of  these  configurations  are  given  by: 

ras  <VS  +  u)a  x  V3)  *  EFg  (9. a) 

and 

[I3]  “s  +  US  X  I  ^S  l“s  *  l-Mg  (9.b) 

where : 

EFg  =  summation  of  the  external  forces  acting  on  the  S/0 
or  S/A 

EMg  =  summation  of  the  external  moments  acting  on  the 
S/0  or  S/A 

ms  =  mass  of  the  S/0  or  S/A 

[Ig]  =  moment  of  inertia  matrix  for  the  S/0  or  S/A 

Vg  =  linear  acceleration  vector  for  the  S/0  or  S/A 

Vg  =  linear  velocity  vector  for  the  S/0  or  S/A 

u>g  =  angular  velocity  vector  for  the  S/0  or  S/A 

cj3  =  angular  acceleration  vector  for  the  S/0  or  S/A 


where: 

3  AC  2 
VAC  = 
®  AC  = 
“AC  = 
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The  poeslble  external  forces  acting  on  the  S/0  are  given  by: 
IF,0  -  Pc  +  ?R  +  Fr  +  PA  +  *D  +  Fdc  +  Frc  +  Fg  ( 10. a) 

where: 

0*0 

Fc  =  catapult  forces 

*0 

FR  =  roclcet  forces 

*0 

Fr  =  rail  forces 

a* 

Fa  =  aerodynamic  forces 

Fd  5  DART  forces 

*0 

Fdc  =  drogue  parachute  forces 

Frc  =  recovery  parachute  forces 

*0 

Fg  =  gravitational  forces 

The  possible  external  moments  acting  on  the  S/O  are: 

EMgo  -  Mc  +  Mr  +  Mr  +  MA  +  MD  +  MDC  +  MRC  ( lO.b) 

where: 


Me 

= 

catapult  moments 

*>0 

Mr 

= 

roclcet  moments 

*0 

Mr 

= 

rail  moments 

*0 

Ma 

= 

aerodynamic  moments 

*0 

Md 

= 

DART  moments 

Mbc 

= 

drogue  parachute  moments 

Mrc 

s 

recovery  parachute  moments 
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The  external  forces  and  moments  acting  on  the  S/A  are: 

M  M* 

zfSA  ■  fA  +  Fg  (11. a) 

EMsa  -  Ma  (ll.b) 


The  new  positions  and  orientations  of  the  S/0  or  S/A  with  respect  to 
the  EFCS  are  obtained  by  integrating  twelve  (12)  differential  equations.  These 
equations  are: 


us 

-  / usdt 

(12. a) 

vs 

•  /vs<it 

(12. b) 

ws 

-  /w3dt 

(12. c) 

xs 

-  Jxsdt 

(12.d) 

ys 

-  /ysdt 

(12. e) 

2s 

-  J*8dt 

(12. f) 

Ps 

■  / PsdC 

(12. g) 

*  /q sdc 

(12. h) 

*s 

-  /rsdt 

(12. i) 

Ps 

-  /psdt 

( 12. j ) 

^s 

•  r4sdt 

(12. k) 

e 

*s 

-  sradc 

(12.1) 

where : 


—  — 

xs 

us 

Fs 

”  ^DSE^T 

vs 

• 

*3 

ws 

— .  _ 

(13) 
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•  •  • 

IC  should  be  noted  that  all  u,  v,  w,  p,  q,  and  r  terms  re£er  to  the  SCS,  while 
the  x,  y,  x,  p,  q  and  r  terms  refer  to  the  E7CS. 

3.3.3  Occupant  Alone.  The  occupant  alone  (0/A)  configuration  cannot  be 

considered  as  a  rigid  body  because  of  the  numerous  degrees  of  freedom  of  the 
various  body  elements.  Consequently,  any  singular  consideration  of  angular 
orientation,  velocity,  or  acceleration  is  meaningless.  Therefore,  the  equation 
of  motion  for  the  0/A  is  based  on  only  the  three  (3)  linear  degrees  of  freedom, 
and  is  given  by: 

“QA^OA  "  zf0A  (14) 

where: 

oqa  =  mass  of  the  0/A 

V0a  S  linear  acceleration  vector  of  the  0/A 
EFqa  =  summation  of  the  external  force  vector  acting  on  the  0/A 
The  external  forces  acting  on  the  0/ A  are : 

IF0A  “  fRC  +  fA  +  Fg  (15) 

The  new  position  of  the  0/A  with  respect  to  the  EFCS  is  obtained  by 
Integrating  the  same  twelve  (12)  equations  as  the  S/0  and  S/A  (equations  12.a-l)  . 
However,  all  angular  terms  are  set  equal  to  zero.  This  causes  the  DCM  to  be  a 
unity  matrix,  and  there  is  no  relevant  distinction  between  the  OACS  and  EFCS 
position  and  velocity  terms.  These  equations  become: 
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U0A  “  J« 0Adt 

(15. a) 

VQA  "  Jv 0Adt 

(15. b) 

WQA  *  JwOAdt 

(15. c) 

X0A  ”  fx OAdt 

(15. d) 

yoA  ■  /yoAdt 

(15. e) 

zOA  "  Sz OAdc 

(15. f) 

POA  ”  0 

(15. g) 

qoA  "  0 

(15. h) 

r0A  “  0 

(15. i) 

POA  "  0 

(15. j) 

^oa  *  0 

(15. k) 

rOA  *  0 

(15.1) 

3.4  Forces  and  Momenta 

The  various  forces  and  moaients  acting  on  the  Individual  system  ele¬ 
ments  result  In  changes  to  the  trajectories  of  these  elements.  These  forces 
and  moments  are  dependent  upon  a  wide  variety  of  independent  variables,  and 
will  be  separately  discussed  in  association  with  their  sources. 

3.4.1  Aircraft.  Prior  to  the  unlocking  of  the  catapult  tubes  at  ejection 
initiation,  the  seat/occupant  (S/0)  is  located  by  a  fixed  position  vector  from 
the  aircraft  center  of  gravity  (CG).  Thus,  the  forces  and  moments  created  by 
the  aircraft's  trajectory  are  directly  translated  to  the  S/0.  However,  since 
these  forces  and  moments  do  not  affect  the  S/0  trajectory  relative  to  the  air¬ 
craft,  the  GESS  program  translates  the  aircraft  accelerations  and  velocities 
directly  to  the  S/0,  but  does  not  calculate  specific  forces  and  moments  exerted 
on  the  S/0  by  the  aircraft.  The  linear  and  angular  aircraft  accelerations  in 
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the  ACS  are  entered  as  input  time  tables  (refer  to  Sections  4. 3. 4.1.  and 
4.3.4. 2.),  allowing  the  user  to  simulate  the  forces  and  moments  of  any  air¬ 
craft  maneuver  (also  refer  to  Section  3.3.1). 

3.4.2  Catapults  and  Catapult  Tube.  The  GESS  program  allows  the  simulation 
of  either  one  or  two  seat  catapults.  Input  variables  allow  the  seat-to-catapult 
moment  arm  vector(s)  ,  rg^,  to  be  specified  according  to  seat  geometry.  If  two 
catapults  are  specified,  the  right  (i  ■  1)  and  left  (i  •  2)  catapult  thrust  time 
histories  are  input  and  tabulated  separately  to  allow  greater  flexibility  in 
simulating  catapult  system  anomalies.  The  input  catapult  thrust  versus  time 
tables  are  linearly  interpolated  to  determine  thrust  magnitudes  between  inter¬ 
mediate  time  points. 


In  the  RCS,  the  right,  left,  and  total  catapult  thrust  vectors  (FCj, 
FC2»  and  Fc,  respectively)  are  given  by: 


fc1  -  co.o, Tcly 

FC2  "  (O.O.Tcz)1 

M  A*  M 

FC  *  FCI  +  fC2 


(16. a) 
(16. b) 


(16. c) 


where  T^  and  T^z  are  the  right  and  left  thrust  table  values,  respectively. 
The  corresponding  catapult  moment  vector.  Mg,  is  formed  to  be: 


^  "  rSCl  x  FC1  ♦  rSC2  *  fC2 


(  '16. d) 


where:  rg^  =  vector  moment  arm  from  catapult  1  line  of  thrust  to  S/0  CG 

Figure  3-6  illustrates  the  forces  and  moments  on  the  seat/occupant  due  to  cata¬ 
pult  thrust. 
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An  additional  force  and  moment  source  sometimes  encountered  in  ejec¬ 
tion  seat  systems  is  that  of  the  bending  catapult  tube.  As  the  seat  moves  out 
of  the  cockpit,  it  will,  in  some  designs,  remain  attached  to  a  single  extended 
catapult  tube  until  the  tube  extension  length  has  been  achieved.  As  the  exten¬ 
sion  of  the  tube  progresses,  the  aerodynamic  forces  acting  on  the  seat/occupant 
cause  the  tube  to  deflect.  This  deflection  results  in  a  countering  force  on  the 
S/0.  Figure  3-7  Illustrates  the  tube  bending  forces  and  moments. 

In  the  RCS,  the  bending  length,  Lb,  is  found  from: 

Lb  -  [(AX)2  +  (AY)2  +  (AZ-Li)2]^2  (17. a) 

where : 

AX  s  change  in  X-axis  position  of  catapult  tube  end 

AY  5  change  in  Y-axis  position  of  catapult  tube  end 

AZ  s  change  in  Z-axis  position  of  catapult  tube  end 

=  initial  catapult  length 


The  catapult  tube  deflection,  6 ,  Is  found  from: 


<5  -  f(AX)2+(AY)2]  1/2 

(17. b) 

The  total  tube  bending  force,  Ft,  exerted  on  the  S/0  is: 

Ft  t 5  (i  -  — )2 

(18) 

where : 

K^.  =  catapult  tube  stiffness  constant 

=  maximum  catapult  tube  extension  length 

This  yields  the  orthogonal  force  components: 

Ftx  -  -FtAX/S 

(19. a) 

Fty  -  -FtAY/5 

(19. b) 

Ftz  "  “UtFt 

(19. c) 

where : 

ut  =  catapult  tube  friction  coefficient 

And,  the  corresponding  tube  bending  moments  are: 
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3.4.3  Sustalner  Rockets.  The  sustainer  rocket  thrusts  are  handled  similarly 
to  the  catapult  thrusts  in  that  each  of  up  to  six  rockets  have  separate  data 
inputs  and  thrust  vs  time  tables,  and  are  computed  separately  within  the  pro¬ 
gram.  The  multiple  inputs  may  be  used  to  investigate  sustainer  rocket  anom¬ 
alies  or  intentionally  introduced  asymmetry  to  achieve  divergence  in  multi-seat 
aircraft.  The  unit  vectors  for  each  rocket's  line  of  thrust,  as  seen  in  the 
SCS,  are  given  by: 

URi  *  (cos  aRi»  cos  ®  Ri »  cos  fRi)*  (i  *  1  to  n,  n  <  6) 


Figure  3-8  illustrates  the  modified  Euler  angles  specifying  the  rocket  thrust 
line  orientation  with  respect  to  the  RCS,  and  the  unit  vector  components  in 
the  SCS. 

The  sustainer  rocket  thrust  vector,  FR,  in  the  SCS  is  given  by: 

Fr  *  ETriUri,  (n  <  6)  (21) 

i»l 

where : 

Tr^  =  rocket  i  total  thrust. 

The  corresponding  rocket  moments  are: 


MR  "i-irR1  *  Trj-UriJ’  (n  *  6)  (22) 

The  effect  of  rocket  fuel  burnoff  on  the  S/0  mass,  rngo.  is  computed 

as  a  linear  function  of  rocket  i  burn  time,  0RBi,  from: 

n  .  0RB4  . 

®S0  *  mS0n  ”**,1  r -  x  mRFi  )  (23) 

0  i-l  OxRBi 

where : 

“SO  =  initial  S/0  mass,  including  rocket  fuel 


®TRBi  =  total  rocket  i  burn  time 


mRFi  =  initial  mass  of  rocket  i  fuel 
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3.4.4  Rail  System.  Several  different  rail  system  configurations  are  used 
for  directing  the  seat/occupant  out  of  the  aircraft.  The  rail  can  be  perma¬ 
nently  mounted  in  the  aircraft  with  either  rollers  or  slider  blocks  (slippers) 
attached  to  the  seat,  or  the  rail  may  be  mounted  on  the  seat  with  slippers 
attached  to  the  aircraft.  There  are  also  ”rail-within-a-rail"  systems  in  which 
both  the  seat  and  the  aircraft  are  provided  rails  which  contact  each  other.  The 
GESS  program  allows  the  user  to  specify  either  type  of  system. 

The  program  determines  the  displacement  of  each  slipper  from  the  rail, 
perpendicular  to  the  RCS  Z-axls ,  at  each  integration  time  step.  For  the  rail- 
within-a-rail  option,  four  slippers,  each  pair  separated  by  the  rail  length, 
are  assumed.  Based  on  this  displacement,  a  restoring  force  is  computed,  con¬ 
straining  the  slider  blocks  to  the  rails. 


Up  to  six  slippers  can  be  specified;  they  are  numbered  consecutively 
in  the  order  in  which  their  initial  positions  are  entered  in  input.  The  de¬ 
flection  vector  from  the  ith  slipper  to  its  respective  rail,  rSSBi,  which  is 
perpendicular  to  the  RCS  Z-axis,  is  denoted  as: 


6SSBi  "  (6SSBix»  5SSBiyt  0)  "  rSSBi  "  6SSBiz  (24. a) 

and  is  illustrated  in  Figure  3-9.  Solving  for  ^gSBi’  we  obtain,  in  the  RCS: 


rSx~ 

rSSBix 

rSBiox 

rSy 

+  DSr 

rSSBiy 

- 

rSBioy 

_  0  _ 

__  0  _ 

_  0  _ 

(24. b) 
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where : 

rSj  = 
DSr  = 
rSSBlj  3 
rSBioj  3 


vector  components  of  S/0  CG  position  in  RCS 
SCS  to  RCS  DCM 

components  of  vector  from  S/0  CG  to  slipper  i  in  SCS 
vector  components  of  slipper  i  initial  position  in  RCS 


The  reaction  forces  at  slipper  i  are  computed  by: 


FSBix  *  “kx  5SSBix 

(25. a) 

FSBiy  *  ~ky  5SSBiy 

(25. b) 

where : 

kj  =  slipper  i  spring  constant  in  j  direction 
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An  additional  pitching  moment,  MT,  caused  by  the  torsional  bending  of  the 
rail  as  the  S/0  angular  position  varies  relative  to  the  RCS,  is  given  by: 

-  -K-r(eSr  -  ®SrO)  <29> 

where : 

iCf  =  torsional  bending  constant 

@gr  =  pitch  angle  of  seat  to  rail 

0grO  =  initial  pitch  angle  of  seat  to  rail 

Therefore,  the  rail/ slipper  moment  vector  acting  on  the  S/0  in  the  SCS,  is: 

6  _  _ 

Mr  -  DSr  l  MSBi  +  Mt  (30) 

i-1 

3.4.5  Aerodynamic  Drag.  The  forces  acting  on  the  seat/occupant  and  seat 
alone  due  to  aerodynamic  drag  are  based  upon  tabular  data  similar  to  that  found 
in  Reference  6.  These  data,  in  the  form  of  drag  coefficients,  are  functions  of 
the  angles  of  attack  and  sideslip  and  of  Mach  number.  The  coefficients  are 
modified,  when  appropriate,  to  compensate  for  the  partial  wind  exposure  encount¬ 
ered  as  the  S/0  begins  to  emerge  from  the  aircraft. 

First,  the  velocity  of  the  S/0  or  S/A  with  respect  to  the  air, 
Usa»  found  from: 


USA  “  US  "  DSE^ 


(31) 
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and  its  scalar  magnitude,  Vg^,  is: 

VSA  *  <USA  +  VSA  +  WSA>1/2 


(32) 


where: 

A# 

dSE 

aw 

me 

USA»  VSA»  WSA 


S/0  or  S/A  velocity  vector  wrt  EFCS 
SCS  to  EFCS  DCM 
wind  velocity  vector 
vector  components  of  UgA 


The  angles  of  attack,  aa,  and  sideslip,  83,  as  shown  in  Figure  3-10,  are 
determined  by: 


°3  *  -can*1  (wSA/uSA)  (33) 

6S  -  sin"1  (vSA/VSA)  (34) 


The  Mach  number  is  found  from: 


where : 

Vx  -  49.0212  x  T1^ 

Tgjjb  =  ambient  air  temperature 

AW 

The  aerodynamic  force  vector,  FA  is  then  given  by: 


(35) 


(36) 


*A  -  PambVsWFs/2 


(37) 
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where: 


Paab  =  ambient  air  density 

fx  =  fraction  of  S/0  exposed  to  the  air  stream 


Ag  =  aerodynamic  reference  area  of  S/ 0  or  S/A 


CFS 


vector  of  S/0  or  S/A  aerodynamic  force  coefficients,  func¬ 
tions  of  as,  Ss,  and 


The  aerodynamic  moment  vector,  M^,  is  given  by: 


M* 


(PambVSArSSCG 


CFS  +  LRSCMs)/2 


(38) 


where : 


rSSCG  =  vector  location  of  the  seat  ceuter  of  gravity  in  SCS 
»>» 

CMS  5  vector  of  S/0  or  S/A  aerodynamic  moment  coefficients,  func¬ 
tions  of  as,0s,  and  Njj 

LRg  s  aerodynamic  reference  length  of  S/0  or  S/A 

As  the  S/0  moves  from  the  cockpit  into  the  airstream,  the  reference 
area,  reference  length,  and  center  of  aerodynamic  forces  are  modified  to  account 
for  the  partial  wind  exposure.  Hence,  if  we  define  4;  as  the  fraction  of  the  S/ 0 
exposed : 


"RSB 


+  ZSA  "  ZCPH 


ZSA 


then: 


4>asf.  (IzrsbI+zsa)  <  ZCP 

asf*  ( IzRSbI+ZSa)  >  ZCP 


(39., 


(39. 
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Similarly: 

’I'I'RSF*  (IZRSBI+ZS)  ^  ZCP  (40. a) 

LrS  “  { 

lRSF*  < |ZrsbI+ZS>  *  zCP  (40. b) 

where: 

Agp  =  aerodynamic  reference  area  of  S/0  fully  out  of  aircraft 

Lrsf  =  aerodynamic  reference  length  of  S/0  fully  out  of  aircraft 

Zrsb  =  travel  of  RCS  Z-position  of  seat  bottom 

ZSA  =  S/A  height  in  RCS 

ZCPH  =  cockpit  height  in  RCS 

The  partially  exposed  aerodynamic  reference  point  (center  of  aerodynamic  forces), 
ZSRP>  is  from: 

ZSRP  “  ZSSRP  +  (ZSBOT  ~  ZSSRP  +  ZSA>  t1  “  *>  (41) 

where: 

ZSSRP  =  fully  exposed  aerodynamic  reference  point 
ZSBOT  5  SCS  Z-position  of  seat  bottom 

As  the  fully-exposed  S/ 0  or  S/ A  rotates,  a  retarding  moment  tends  to 
oppose  this  rotation  due  to  aerodynamic  damping.  GESS  has  incorporated  this 
moment  in  the  form  of  a  damping  factor,  fp,  such  that: 

a#  a# 

“Sd  "  “S  fD  (42) 


where : 
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and 

(Dg^  =  dampened  rotational  velocity  of  S/0  or  S/A 

*s# 

uig  =  calculated  rotational  velocity  of  S/0  or  S/A 
Cj)  =  input  damping  coefficient  (0<C[j<l) 

At^  =  integration  time  step  for  simulation  phase  i 


In  the  OACS,  the  aerodynamic  forces  are  computed  based  upon  a  fixed, 
multidimensional,  aerodynamic  reference  area  where  no  rotational  effects  are 
assumed.  First,  the  0/A  velocity  with  respect  to  the  air,  UQA,  is  found  from: 


uOA  M  u0  "  WE 
and  its  scalar  magnitude,  VqA,  is: 

2  2  2  II? 

VOA  3  <U0A  +  V0A  +  wOA^ 


(44) 


(45) 


whe  re : 


Uq  =  0/A  velocity  vector  wrt  EFCS 

U0A>  v0A>  w0A  2  vector  components  of  Uqa 


The  O/A  angles  of  attack,  a0A,  and  sideslip,  8qA,  are  f°und  from: 

a0A  *  sin"l(w0A/V0A)  (4,,) 

80A  “  tan_1^v0A/u0A)  (i7> 

The  vector  components  of  the  aerodynamic  force  coefficient,  Cp,  are  determined 
from: 


CFX  3  -SIGN[(cosa0Acos60A)»  u0A] 

(48. a) 

CFY  *  -SIGM[ (cosa0AsinS  0A) ,  v0A) 

(48. b) 

CFZ  -  -SIGN[(sina0A),  w0A] 

( 48  .  c  ) 
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where  Che  SIGN  function  assigns  Che  firsc  argument  Che  sign  of  Che  second. 
Therefore,  Che  Cocal  veccor  OA  aerodynamic  force,  FA0,  is 

”*AO  “  PambVOAAQA(V2  <49) 

3.4.6  Stabilization.  Various  mechanisms  are  used  to  stabilize  the  S/O 
during  sustainer  rocket  burn.  The  GESS  program  includes  provisions  for  simu¬ 
lating  a  Directional  Automatic  Realignment  of  Trajectory  (DART)  system.  As 
the  seat  moves  away  from  the  aircraft,  the  DART  lines  play  out  through  a  brake 
mechanism  in  the  DART  frame.  The  brake  force  is  assumed  to  act  along  each 
DART  line  from  the  time  that  the  DART  line  first  becomes  stretched  to  the  time 
that  the  lines  separate  from  the  seat  (see  Figure  3-11). 


The  vector  DART  force,  FDi,  for  a  given  DART  line  i,  presented  in  the 
SCS,  is  given  by: 


where 


rDi 


1  kdf  ~ 

-  rDCAPi 

2  LDLt 


(50) 


kDF 


DART  force  constant 


length  of  dart  line  i  from  confluence  to  attachment 
points 


rDCA?1 


vector  from  confluence  to  attachment  points  of  dart 
line  i . 


The  total  DART  force,  Fq,  then  becomes: 


FD 


(51) 


The  vector  DART  moment,  MD1 ,  for  a  given  DART  line  i,  in  the  SCS,  is 


given  by: 

MDi  "  rSCPi  x  FDi 


(52) 
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where: 


vector  from  the  seat  center  of  gravity  to  the  aircraft 
confluence  point  for  DART  line  i 


The  total  DART  moment,  Mq,  then  becomes: 


f  «Di 

i-1 


(53) 


3.4.7  Drogue  Slug/Container.  The  aerodynamic  drag  force  applied  to  the 

drogue  slug/container ,  ?DSC»  after  its  deployment  is  used  to  determine  the 
drogue  parachute  acceleration.  This  force  magnitude  is  determined  by: 

FDSC  “  PambCDSCADSCVDA/2  (54> 

where : 

Cose  =  drogue  slug/container  drag  coefficient 

Aqsc  =  drogue  slug/container  effective  drag  area 

Vqa  =  velocity  magnitude  of  the  drogue  parachute  being  deployed 

by  the  slug/container  relative  to  the  air 

2  2  2  1/2 
"  <UDA  +  VDA  +  WDA> 

UDA*  VDA»  WDA  =  components  of  vector  Uqa 

Uda  =  dseUd“W 

Dse  =  SCS  to  EFCS  DCM 

Uq  =  drogue  parachute  velocity  vector  in  EFCS 


W 


=  wind  vector  in  EFCS 
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The  axial  components  of  Fpgg  are  found  from: 


tp  m  DSC  rr 

*DSC  -  °DA 

VDA 


(55. a) 


and  Che  initial  slug/conCainer  acceleration,  ,  becomes; 

o 


aDSCo“  —  FDSC 
“DSC 


(55. b) 


where: 


mQsc  =  mass  of  the  drogue  slug/container 


3.4.8  Parachutes .  Both  drogue  and  recovery  parachute  forces  and  moments  are 
calculated  similarly.  Parachute  release  or  deployment  occurs  at  the  time  the 
parachute  begins  to  leave  the  seat.  After  deployment,  the  parachute  begins  to 
develop  or  play  out,  but  will  not  exert  any  force  until  the  suspension  line  is 
fully  extended  (line  stretch).  Upon  achieving  line  stretch,  the  parachute 
begins  to  fill  and  exert  force.  The  total  time  required  to  completely  fill 
the  parachute  is  nearly  always  a  function  of  the  velocity  of  the  S/0,  and  is 
determined  by  linear  extrapolation  from  a  user  provided  data  table.  The  par¬ 
tially  Inflated  diameter  of  the  parachute,  dn,  as  a  function  of  time,  t,  is 
calculated  as: 


t-t 


‘full 


Is 


Atft 


<cls<  c  <  tls^ft’ 


(56) 


where : 

dfull  =  fully  inflated  parachute  projected  diameter 
ti3  =  time  at  which  parachute  line  stretch  is  achieved 

AtfC  =  total  time  required  for  parachute  fill 
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The  total  force  magnitude  exerted  by  the  parachute  on  the  S/0  or  0/A 
is  found  from: 

Fe  -  -<’a«Pvl*cop4/8  <”> 

where 

VgA  =  total  resultant  velocity  of  S/0  wrt  air 

2  2  2  1/2 
”  ^sa^sa^sa) 

Cqp  2  parachute  drag  coefficient 


parachute 

force 

vector,  Fp,  is  componentized  as: 

Fpx 

-  Fp 

cosap  cos8p 

(58. a) 

FPy 

-  Fp 

cosap  sinSp 

(58. b) 

FPz 

»  Fp 

sinOp 

(58. c) 

where : 

atp  =  parachute  angle  of  attack  wrt  air 

-  sin'1  («SA/VSA) 

8p  2  parachute  angle  of  sideslip  wrt  air 

-  tan'1  <vSA/uSA) 


rgp  2  position  vector  from  S/0  center  of  gravity  to  parachute 
confluence  point 
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3.4.8. 1.  The  drogue  parachuCe  (drogue)  is  deployed  from  Che  seat  by  che  drogue 
slug/coaCalaer  until  che  drogue  suspension  line  is  completely  extended  (line 
stretch).  The  total  time  required  to  achieve  line  stretch  may  be  specified  on 
input  as  either  extrapolation  table  values  based  on  S/ 0  linear  velocity,  or  it 
may  be  calculated  based  upon  the  acceleration  provided  by  the  slug/container. 
At  this  time  the  drogue  begins  to  fill.  The  total  time  required  to  completely 
fill  the  drogue  is  also  determined  by  linearly  extrapolating  input  table  values 
based  on  S/0  velocity. 

An  option  is  provided  to  allow  the  simulation  of  a  duplex  VELCON 
drogue  system.  The  VELCON  drogue  consists  of  two  tandem  drogues,  a  relatively 
small  one  preceded  by  a  larger  one.  The  two  drogues  are  connected  by  a  vel¬ 
ocity-calibrated  breakaway  bridle.  At  high  speeds,  the  force  generated  by  the 
larj;e  drogue  causes  the  bridle  to  fail,  separating  it  from  the  system  and 
limiting  the  force  imposed  on  the  S/ 0.  At  low  speeds,  the  forces  of  the  two 
drogues  are  additive.  The  speed  at  which  the  large  drogue  separates  is  part 
of  Che  program  input,  along  with  the  descriptions  of  che  two  drogues. 

3. 4. 8. 2.  The  recovery  parachute  (chute)  is  deployed  from  the  seat  by  various 
mechanisms,  depending  upon  such  variables  as  Che  ejection  altitude,  velocity, 
drogue  force,  and  others.  The  GESS  program  allows  che  user  to  specify  either  a 
fixed  time  of  deployment  or  to  use  an  extrapolated  time  of  deployment  based  on 
table  values  that  are  functions  of  the  S/0  velocity.  The  user  can  also  specify 
*ne  altitude  above  which  the  chute  will  not  deploy  (CHALT2),  and  the  altitude 
below  which  the  parachute  will  deploy  at  the  specified  deployment  time  (CHALT1). 
Between  these  two  altitudes  the  user  can  specif iy  an  additional  fixed  time  delay, 
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or  a  variable  delay  based  upon  Che  seat  vertical  deceleration.  The  latter  option 
determines  a  time  of  deployment,  t^ep^0y,  based  on: 

tdeploy  "  cdelay  +  (60) 

fixed  deployment  time  delay 

number  of  time  steps  during  which  the  vertical 
acceleration  exceeds  a  specified  acceleration  limit 

integration  step  time  increment 

The  total  time  required  to  achieve  line  stretch  is  specified  in  input 
as  linear  extrapolation  table  values,  based  upon  S/0  linear  velocity.  The  user 
has  the  option  of  determining  parachute  fill  time  by  extrapolation  tables  based 
on  S/0  velocity  or  by  calculations  based  on  theory.  The  latter  approach  calcu¬ 
lates  the  area  under  the  curve  representing  current  parachute  volume,  and  then 
compares  this  volume  to  the  actual  parachute  volume  at  full  inflation.  If  the 
two  volumes  vary  by  more  than  a  given  tolerance  limit,  a  new  value  for  parachute 
filling  time  is  assumed,  and  the  area  under  the  new  volume  curve  is  calculated. 
The  proper  value  for  parachute  filling  time  is  determined  when  the  two  volumes 
are  within  the  tolerance  limit.  A  more  detailed  discussion  can  be  found  in 
Reference  7. 

The  maximum  parachute  volume,  ^max  *  is  found  from: 


where : 


cdelay 

i 


At 


where: 


(61) 


Dp 


projected  parachute  diameter 
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The  projected  parachute  area,  Ap,  is  found  from 


AP  “  ?APi 
1-1 


where : 

ZiRu 


3(P2R2i+P3i/A) 


Zi 


TDi 

TDi 

td± 


.9 


Rli  "  TDi(PllTDl  _2*2p) 
r21  “  P31  lo8e  (p3i/A)  “  11»25 

pli  "  1  +  TD1(2.2p-l) 

r21  “  BVq/253. 125 

-  11.25  Tqj[  +  A 

Tut  -  fraction  of  total  chute  filling  time 

-  0.  to  1.0 


5000  W0A 

g(p/p0)Dooc 


1 20TppCjjSp 
B  -  — j - 

dodc 

Wqa  =  weight  of  the  occupant  alone  after  seat  separation 

g  =  gravitational  acceleration  constant 

p  s  localized  air  density 

pa  2  air  density  at  sea  level 
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°ODC 

Dp 

TpP 

CD 

SP 

P 

*o 


parachute  flattened  diaater  -  */ Dp 
parachute  projected  diameter 

estimated  parachute  filling  times.  Initially  -  1.0  sec 

parachute  drag  coefficient 

o 

parachute  projected  area  -  tt Dp/ 4 

parachute  porosity 

free  stream  air  velocity 


The  calculated  parachute  volume  becomes : 

2 

tFP^od0DCaC 

Vol  -  -  (64) 

w 


This  volume  is  then  compared  with  the  actual  volume: 


0.01 


(65) 


If  the  tolerance  condition  is  not  met,  then  Tpp  is  incremented  or  decremented  by 
half  of  the  difference  between  the  curren  value  and  either  its  closest  previous¬ 
ly  tried  value  or  zero. 

3.5  Special  Features 


3.5.1  Dynamic  Response  Index.  The  user  is  offered  the  option  of  having  the 
program  determine  the  maximum  value  of  the  Dynamic  Response  Index  (DRI)  during 
the  time  period  of  seat  motion  prior  to  rail  separation.  The  DRI  is  a  standard 
measure  of  the  acceleration  stress  applied  to  the  seat  occupant  during  ejection. 
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The  DRI  Is  determined  at  every  integration  time  step  during  in  air¬ 
craft  acceleration  by  first  numerically  integrating  the  following  second  order 
equation  given  in  Section  6.5.1  of  MIL-S-9479B(USAF) : 


dSz  -  VZ 

dVz  -  Az  -  23.7  Vz  -  2798. 41SZ 


az  -  z3  ~  (z.  -  Zs.iJUt-t.pMt) 


where: 


AZ 


“s 

Vi  5 

At  i 


vertical  acceleration  in  SCS  during  time  step 
vertical  acceleration  in  SCS  at  time  t 
vertical  acceleration  in  SCS  at  time  t_^ 
length  of  integration  time  step 


Hence,  the  time  step  value  of  DRI,  Iqr,  is  taken  from: 


XDR  “  CDRISZ 


where : 


Cqrx  =  dynamic  response  index  constant 


and 


Sz  -  JdSz 

-  IJdVz 


(66. a) 
(66. b) 


(66. c) 


(67. a) 


(67. b) 
(67. c) 


When  selected,  the  DRI  option  will  save  the  maximum  DRI  value  calcu¬ 
lated  over  the  pre— separation  period  and  the  time  at  which  this  value  occurred. 


3 


•vj 

% 

\ 
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3.5.2  Dynamic  Center  of  Gravity.  Previous  ejection  seat  tests  have  demon¬ 
strated  a  considerable  oscillatory  motion  between  the  occupant  and  the  seat 
prior  to  seat/occupant  separation.  Since  the  relationship  between  the  S/0 
center  of  mass  and  the  rocket  thrust  line  can  be  critical  to  the  accurate 
determination  of  the  S/0  trajectory,  a  user-optional  routine  has  been  included 
in  CESS  to  approximate  the  shifting  of  the  S/0  center  of  gravity  (CG)  location 
throughout  this  phase  of  the  simulation.  The  development,  methodology,  and 

evaluation  of  this  dynamic  CG  mathematical  model  used  for  this  routine  is 
described  in  Reference  8. 

The  dynamic  CG  model  is  represented  by  an  orthogonal  three  degrees  of 
freedom  spring  and  damper  system.  A  two-dimensional  representation  of  this 
system  is  shown  in  Figure  3-12.  The  general  equation  of  motion  for  such  a 
system  is  given  by: 

M  M  IW 

5+  2cui0d  +  wjs  -  A(t)  (68) 

where 

6  =  deflection  vector  of  occupant  relative  to  seat 

•  ** 

6,5  =  first  and  second  time  derivatives  of  deflection  vector 

c  - 

2mW0 

b  =  damping  constant 

m  =  mass  of  occupant 

u0  »  (K/m)l/2 

K  =  spring  constant 

A(t)  =  time  dependent  external  acceleration 

A  major  feature  of  the  dynamic  CG  model  is  its  use  of  non-linear  spring 
constants.  The  spring  constants  of  all  three  axes  incorporate  different  values 
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be c ween  Che  positive  and  negative  displacement  directions.  Both  the  SCS  X  and  Z 
axes  incorporate  a  positive  "slack"  displacement  zone  in  which  there  is  no  spring 
force  exerted  on  the  occupant.  Further,  the  Z-axis  incorporates  two  separate 
spring  constant  values  in  the  negative  direction;  this  simulates  the  effects  of 
the  initial  compression  of  the  seat  cushion  before  reaching  the  "bottoming"  dis¬ 
placement  corresponding  to  the  cushion  thickness.  Hence,  the  vertical  (Z-axis) 
spring  constant  can  be  represented  as : 

Kg  *  0,  Z  ■  0  to  Zsiack 

H  -  V  ^  ^slack  to  +*“ 

KZ  “  Z  ■  0  to  Z|Jottom 

KZ  “  Kn2  >  Z  ”  Zj,ottom  to 

Typical  dynamic  CG  spring  constant  curves  are  presented  in  Figures  4-16,  4-17, 
and  4-18  on  pages  4-35  to  4-37. 

The  displacement  of  the  occupant  CG  relative  to  the  seat  is  determined 
by  numerically  integrating  Equation  68  at  each  time  step  until  seat/occupant  sep¬ 
aration.  This  displacement  is  then  used  to  update  the  S/0  CG  location  in  the 
SCS,  upon  which  the  remaining  position  vectors  are  similarly  updated. 

3.5.3  Thrust  Vector  Control.  The  user-optional  Thrust  Vector  Control  (TVC) 
feature  simulates  the  vertical- seeking  maneuver  Incorporated  into  the  develop¬ 
mental  Maximum  Performance  Ejection  Seat  (MPES).  In  this  manuever,  a  single 
seat  booster  rocket  (SBR),  located  on  the  seat  bottom,  can  be  glmballed  to 
control  the  roll  and  pitch  movement  of  the  S/0  in  a  manner  which  attempts  to 
bring  the  S/0  to  an  upright  orientation.  An  on-board  microprocessor  updates 
the  direction  cosine  matrix  (DCM)  from  sensor  data  identifying  the  current 
seat  orientation  at  fixed  time  Intervals  from  catapult  ignition  to  rocket 
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burnout.  A  series  of  TVC  thrust  signal  transfer  equations  are  numerically 
Integrated  to  produce  the  orientation  Instruction  given  to  the  rocket  nozzle 
rotation  actuators.  The  development  of  the  vertical-seeking  maneuver  algorithm 
upon  which  the  TVC  routine  is  based  was  conducted  at  the  Naval  Weapons  Center , 
China  Lake,  CA.  The  brief  formulation  that  follows  has  been  included  for  com- 
pleteness  only.  A  more  detailed  analysis  than  that  which  follows  can  be  found 
in  Reference  9. 

The  equations  governing  the  TVC  pitch  rate,  compensation  are  as  follows: 


*3 

«  5OO(0ST-X3) 

(69) 

X4 

-  x5 

(70) 

x5 

-  250000 (X3-X4) -5OOX5 

(71) 

X© 

-  C2952  -  C30  X4  -100X6 

(72) 

-  5(X6  +  20  X6),  |©x|  Olirnit 

(73) 

where: 


= 

signal  transfer  variables 

Xi 

= 

time  derivative  of  X^ 

• 

©ST 

= 

pitch  rate  of  the  TVC  coordinate  system  (TVCCS) 
wrt  the  SCS 

5i 

= 

rocket  command  angles 

Ci 

= 

fixed  constants 

©T 

2 

intermediate  value  of  rocket  pitch  wrt  TVCCS 

* 


rocket  pitch  angle  limit 


©limit 
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The  corresponding  TVC  roll  compensation  equations  are: 


*7 

-  500(*SX-X7) 

(74) 

*8 

-  X9 

(75) 

*9 

-  250000(X7-X8) -5OOX9 

(76) 

*10 

*  C31C1-C32X8-IOOX10 

(77) 

*T 

“  5(Xio+20X1Q) 1  I^Tl^^limit 

(78) 

where 

e 

*ST 

= 

roll  rate  of  the  TVCCS  wrt  the  SCS 

♦t 

= 

intermediate  value  of  rocket  roll  wrt  TVCCS 

♦limit 

= 

rocket  roll  angle  limit 

The  equations  governing  control  of  the  rocket  pitch  angle  are  then: 

*11 

-  2861000  Ox-31. 45XU 

(79) 

e 

0.0,  xx  <  <ellmlt  +  0.25) 

*1 

■ 

(80) 

_x12*  <0Umit  +  0.25)  <  X! 

*12 

-  XU  -  48.26X12  -  90992XX 

(81) 

«1 

-  Xi  +  eRT 

(82) 

and  the  corresponding  rocket  roll  angle  control  equations  are: 

*13 

-  28610009x-31.45X13 

(83) 

*2 

1  0.0,  X2  <  (*liait  +  0.25) 

(84) 

Lx14»  (♦limit  +  0.25)  <  X2 

• 

x14 

-  Xl3-48. 26X^-90992X2 

(85) 

«2 

-  X2 

(86) 
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where: 

Qi  =  rocket  thrust  line  angles,  generated  by  the  signal 
transfer  equations,  vrt  TVCCS 

0gT  =  initial  angle  of  rocket  wrt  TVCCS 

We  now  find  the  TVC  Euler  angles  to  be: 

°TVC  m  s*n  ^1  cos  ^2 

0TVC  “  s^n 

YTVC  “  cos  ^1  cos  ^2 


(87. a) 
(87. b) 
(87  .c) 


Next,  the  TVC  microprocessor  code  is  simulated  by  integrating  the 
angular  rate  data  as  follows: 


di,j 

-  Xidtj  (i-1,2,3;  j-1,2) 

(88) 

Tdi 

-  1  dt  j 

j-1 

(89) 

AC 

“  Cs5t 

(90) 

If  |Tdil 

>  Ac, 

then  d1>2-  .007812Tdt/|Tdi| 

(91.a) 

and  difl-  Tdl-.0078125Tdi/|Tdi| 

(91 .b) 

If  |Tdil 

<  AC, 

then  dlf2 -  Tdi 

(92. a) 

and  ddj^»  0 

(92. b) 

where : 


dI»  j 

= 

angular  displacement  matrix 

= 

vector  time  derivative  of 

TVC 

angular  velocity  wrt  SCS 

«tj 

= 

Integration  time  step  length 

for 

increment  j 

Tdi 

= 

as  defined 

AC 

= 

as  defined 

Cs 

TVC  processor  sampling  rate 
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Using  Che  Crowder-Hession  algorithm,  the  direction  cosine  matrices 
are  updated  as  follows: 

<*2  “  “  ^3,2  ^1  (93. a) 

*^3  *  ^3  "**  ^2,2  ^1  (93. b) 

dl  *  di  +  dj^  ^2  (93. c) 

^3  *  ^3  ”  ^1,2  ^2  (93. d) 

dl  “  d(  —  ^2,2  ^3  (93. e) 

*  ^2  +  ^1,2  ^3  (93. f) 

Finally,  the  rocket  commands  are  updated,  depending  upon  the  DCM 
values,  as  follows: 

If  d3>0: 

51  •  ~d2  (94. a) 

52  *  dl  (94. b) 

If  d3<0,  dj^  ■  d2  *  0,  and  0gx  *  0: 

51  -  2  (95. a) 

52  ■  0  (95. b) 

If  d3<0,  dj  ■  d2  ■  0,  and  0gj  j*  0: 

5 i  “  2  0ST/|0gT|  (96. a) 

52*0  (96. b) 


If  d3<0,  and  ! d  ^ | < | d 2 1 1 
5l  -(  I d 2 1  ”2)d2/  i d 2 1 


(97. a) 
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If  d3<0,  and|d1|>|d2| : 

51  -  -d2 

52  -  (2-ldjl)  <i  1/  | d i | 


(98. a) 
(98. b) 


3.6  Numerical  Integration 

Solution  of  the  various  differential  equations  found  in  all  sections 
of  the  GESS  program  requires  numerical  integration  techniques.  Two  different 
types  of  integration  methods  are  used:  the  Runge-Kutta  method,  and  the  predic¬ 
tor-corrector  method*  The  Runge-Kutta  method  is  well  known  for  its  accuracy, 
but  is  not  computationally  efficient  for  programs  with  numerous  time  steps  and 

large  arrays.  The  predictor-corrector  method  yields  reasonable  accuracy  with 

% 

excellent  computational  efficiency. 

3.6.1  Runge-Kutta.  Runge-Kutta  formulas  involve  the  computation  of  the 
value  of  a  multivariate  function  at  specified  points  without  the  prior  calcu¬ 
lation  of  the  higher  derivatives  of  the  function.  If  we  assume  such  functions 
as: 


-  f(t,x,y) 

(99. a) 

dt 

iZ.  -  g(t,x,y) 

(99. b) 

dt 
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we  will  seek  solutions: 

x(t«i)  -  x(t)  +  Ax 

and 

y( c+h)  -  y(t)  +  Ay 

where: 

h  5  time  step  Interval  »  At 
Ax  *  aki  +  bk2  +  ckj  +  dk4 

Ay  ”  a*l+b*2  +  cJl3  +  <^4 

and  the  general  form  of  the  solution  is  found  from: 


(100. a) 


( 100. b) 


(101. a) 
(101. b) 


ki  -  hf (t,x,y) 

(102. a) 

k2  “  hf(t  +  mh,  x  +  mk j ,  y  +  mij) 

(102. b) 

k3  ■  hf[t  +  nh,  x  +  nk2  +  (n  -  r)ki. 

y  +  ni2  +  (n  -  r)ijJ 

( 102. c) 

k4  ■  hflt  +  ph,  x  +  sk2  +  tkj  +  (p-s-t)ki. 

y  +  si  2  +  tJt3  +  (p-s-tlii) 

( 102. d) 

and  similarly; 


*1 

-  hg(t,x,y) 

(102. e) 

i2 

■  hg(t  +  mh,  x  +  mkj,  y  +  mij) 

(102. f) 

*3 

■  hg[t  +  nh,  x  +  nk2  +  (n-rjkj,  y  +  ni2  +  (n-r)i]J 

(102. g) 

*4 

•  hg[t  +  ph,  x  +  sk2  +  tk3  +  (p-s-t)ki. 

y  +  si  2  +  ti3  +  (p-s-tlij 

(102. h) 

The  method  may  be  extended  for  as  many  first-degree  differential  equations  as  are 
required.  A  more  detailed  discussion  of  Runge-Kutta  methods  can  be  found  in  any 


numerical  methods  textbook 
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There  are  several  variations  of  the  Runge-Kutta  method,  each  with  Its 
own  advantages  and  disadvantages.  In  GESS,  two  variations  are  used.  The  first, 
the  Kutta-Simpson  one-third  rule  variation,  is  used  for  Integrations  involving 
relatively  small  arrays  and  short  simulated  time  periods,  such  as  In  the  calcu¬ 
lations  of  the  DRI  (Section  3.5.1)  and  dynamic  CG  (Section  3.5.2).  This  method 
assumes  constant  values  of : 


m 

■  n  ■ 

r  -  1/2 

(103. a) 

P 

-  t  - 

1 

(103. b) 

s 

-  0 

(103. c) 

a 

-  d  - 

1/6 

( 103. d) 

b 

-  c  * 

1/3 

(103. e) 

GESS  uses  the  Gill  variation  for  integrating  the  very  large  trajec¬ 
tory  arrays  during  the  initial  four  time  steps,  since  it  minimizes  the  number 
of  memory  locations.  This  method  assumes  constant  values  of: 


m 

m 

n  -  1/2 

(104. a) 

P 

m 

l 

( 104. b) 

r 

m 

1//2 

(104. c) 

s 

m 

-1//2 

(104. d) 

t 

m 

1  +  1//2 

(104. e) 

a 

m 

d  -  1/6 

(104. f) 

b 

m 

(1-1//2  )/3 

(104. g) 

c 

m 

(1+1//2  )/3 

(104. h) 
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3.6.2  Predictor-Corrector .  Predictor-Corrector  forward  integration  methods 
involve  the  use  of  an  initial  prediction  formula  for  the  next  time  step  value  of 
a  function,  followed  by  the  application  of  a  more  accurate  corrector  formula.  By 
iterating  these  formulas  until  the  solution  converges,  this  procedure  achieves 
successive  improvements  in  both  the  function  value  and  its  subsequently  calcu¬ 
lated  derivative.  However,  the  Hamming  variation  of  applying  the  predictor- 
correction  method  eliminates  the  need  for  iteration  by  further  modifying  each 
of  the  two  equations  based  on  a  previously-derived  error  estimate ( *0) .  Thus,  the 
two-pass  equations  used  in  GCSS  based  on  previously-determined  values  of  the 
function  and  derivative  at  equally  spaced  time  intervals,  become: 

(105. a) 
(105. b) 
(105. c) 
(105. d) 

where: 

yit  =  value  of  function  at  time  step  k 

y£  =  value  of  function  derivative  at  time  step  k 

p^  =  predicted  value  of  function  at  time  step  k 

cjt  =  corrected  value  of  function  at  time  step  k 

h  s  time  step  interval 


•  Predictor:  p+1  -  y_3  +  ith(2y{j  -  ylL  +  2y'2) 

•  Modifier:  y  -  p+1  -  iH<p0-c0) 

p  121 

•  Corrector:  c+1  -  i(9yg  -  y_2  +  3h(y4.L  +  2yg  -  ylL) ] 

8 

e  Final  Value:  y+1  -  c+1  +  -2-  (p+1  -  c+j,) 


.•.V.VlslVAV'.s .  • 
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3.6.3  Integration  Array.  The  large  number  of  equations  and  saved  values 
associated  with  many  of  the  trajectory  integrations  in  GESS  resulted  in  the 
development  of  a  highly  structured  integration  array  logic.  This  logic  main¬ 
tains  the  current  value  of  each  integration  equation,  in  addition  to  all  inter¬ 
mediate  or  previous  equation  values  for  both  the  Runge-Kutta  and  predictor- 
corrector  methods.  Thus,  the  initial  four  time  steps  can  be  integrated  by  the 
Runge-Kutta  routine,  and  subsequent  integrations  performed  by  the  predictor- 
corrector  routine. 

The  first  value  of  each  array  indicates  the  number  of  equations  in¬ 
volved  in  the  integration.  The  remaining  array  values  are  then  ordered,  in 
accordance  with  the  number  of  equations,  as  described  in  Table  3-1. 
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Array  Location 
1 

2  —  <  1+N) 

(2+N)  —  (1+2N) 

(2+2N)  —  (1+3N) 
(2+3N)  —  (1+4N) 
(2+4N)  —  ( 1+5N) 
(2+5N)  —  (1+6N) 
(2+6N)  —  ( 1+7N) 
(2+7N)  —  (1+8N) 
(2+8N)  —  (1+9N) 
(2+9N)  — (1+10N) 
(2+10N)—<  1+1  IN) 
(2+llN)—(  1+12N) 
<2+12N)—< 1+13N) 
(2+13N)— (  1+14N) 
(2+14N) — ( 1+15N) 
(2+15N)—(  1+1 6N) 


Table  3-1.  Integration  Array  Values 

Description 
number  of  equations  (N) 
results  of  integrations  -  equation  values 
equation  derivatives 
current  Runge-Kutta  K-values 
intermediate  Runge-Kutta  K-value  summations 
intermediate  equation  values 
equation  values  at  (tlme-At) 
equation  values  at  (time-2At) 
equation  values  at  (time-3At) 
equation  values  at  (time-4At) 
equation  derivatives  at  (tlme-At) 
equation  derivatives  at  (time-2At) 
equation  derivatives  at  (time-3At) 
predictor  values  at  currect  time 
predictor  values  at  (time+At) 
corrector  values  at  current  time 


corrector  values  at  (time+At) 
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4.0  INPUT  PREPARATION 


4.1  General  Procedures 

The  GESS  program  has  been  developed  Co  use  a  series  of  specific  steps 
or  procedures  for  its  proper  execution.  The  program  was  developed  using  the 
Control  Data  Corporation  computers  located  at  the  Naval  Air  Development  Center 
with  the  KRONOS  operating  system,  and  was  programmed  using  the  ASCII  FORTRAN  IV 
language^!*  .  Slight  modifications  may  be  necessary  if  the  program  is  to 
be  used  on  other  computer  systems. 

The  first  step  in  the  execution  of  the  program  is  to  create  a  series 
of  aerodynamic  coefficient  tables  using  a  pre-processing  program  called  ACT(l^). 
This  routine  analyzes  a  file  of  raw  data  from  wind  tunnel  tests  consisting  of 
aerodynamic  forces  and  moments  acting  on  the  seat /occupant ,  occupant  alone,  or 
seat  alone.  The  program  generates  a  random  access  file  of  numerically  smoothed, 
tabular,  aerodynamic  coefficients  which  are  used  in  subsequent  runs  of  the  GESS 
program.  Once  a  table  has  been  generated  for  a  specific  seat  or  occupant 
configuration,  it  may  be  reused  for  all  subsequent  GESS  runs  relevant  to  that 
configuration.  The  procedure  for  generating  the  ACT  file  is  generally  discussed 
in  Section  4.2  and  detailed  in  Appendix  A.  A  description  of  the  theoretical 
concepts  involved  can  be  found  in  Reference  14. 

Next,  the  GESS  input  file  is  created.  All  program  variables  are 
listed  in  a  specific  format,  which  is  detailed  in  Section  4.3.  If  a  previous 
input  file  exists  for  the  specific  seat  and  occupant  configuration,  a  relatively 
few  variable  changes  can  usually  produce  the  new  input  file  with  minimum  effort. 
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Finally,  Che  simulation  job  scream  file  is  creaCed.  This  is  Che  lisc 
of  commands  Co  Che  compuCer  which  will  concrol  Che  program  execuCion  and  assem¬ 
ble  Che  ouCpuC  report  daCa.  The  job  scream  is  decailed  for  the  KRONOS  operating 
sysCem  in  Section  4.4. 

After  Che  submission  of  Che  job  scream  and  subsequent  program  execu¬ 
Cion,  a  series  of  simulacion  reports,  detailed  in  Section  5.1,  are  available 
from  Che  outpuC  file.  These  reports  may  be  reviewed  and  analyzed  according  to 
Che  user's  specific  needs.  Also,  two  plotting  files  may  be  created  at  Che 
user's  option  for  further  analysis  using  post-processing  routines,  such  as  Che 
DRAS  utility  program.  This  program,  discussed  in  Section  5.2,  and  detailed  in 
Reference  15,  performs  various  data  manipulation  and  graphic  function  routines 
using  the  data  from  the  program  plotting  files.  The  CALCOMP-type  plots  gener¬ 
ated  can  greatly  enhance  the  interpretation  of  program  output . 

4.2  Aerodynamic  Coefficient  Tables 

During  the  execution  of  the  GESS  program,  equally  spaced  aerodynamic 
coefficients  are  required  for  computing  the  aerodynamic  forces  acting  upon  Che 
seat/occupant  system.  These  coefficients  are  usually  derived  through  wind 
tunnel  experiments  using  the  subject  ejection  seat(^).  The  forces  and  moments 
data  produced  by  these  tests  are  used  as  input  for  the  Aerodynamic  Coefficient 
Table  (ACT)  program(^),  which  creates  the  necessary  coefficient  tables  as 
functions  of  the  seat  orientation  and  velocity.  The  ACT  program  provides  an 
efficient  method  for  creating,  modifying,  and  storing  these  aerodynamic  coeff¬ 
icient  tables  on  random  access  files  for  retrieval  by  GESS  during  simulation 
execution.  The  ACT  program  was  closely  modeled  after  the  RFWTHR  program  (^). 


v:-.‘ 
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Detailed  instructions  for  using  the  ACT  program  are  given  in  Appendix 
A.  These  instructions  have  been  extracted  from  Reference  14  because  of  the  impor¬ 
tance  of  having  proper  aerodynamic  coefficient  tables  for  the  successful  execu¬ 
tion  of  GESS. 

4.3  Simulation  Input  Variables 

Prior  to  running  the  GESS  program,  an  80  column  card  image  input  file 
must  be  created.  This  file  consists  of  START  and  STOP  cards,  a  three  card  head¬ 
ers,  and  fifteen  data  sections.  A  sample  listing  of  a  typical  input  file  is 
given  in  Appendix  B. 

The  program  can  be  run  against  several  sets  of  data  without  the  user 
having  to  resubmit  the  job.  This  is  accomplished  by  placing  the  START  card  at 
the  beginning  of  each  input  sec  so  that  the  input  sets  are  separated  by  START 
cards.  A  STOP  or  blank  card  is  placed  at  the  end  of  the  entire  input  file.  The 
job  is  terminated  when  the  program  reaches  this  instruction. 

The  three  line  header  is  used  to  identify  the  input  file,  as  well  as 
all  subsequent  output  reports.  Any  alphanumerical  data  may  be  entered,  but  the 
third  line  serves  to  identify  the  seat  to  be  simulated.  If  three  lines  (or 
records)  of  description  are  not  necessary,  then  blank  records  must  be  inserted 
to  complete  the  mandatory  three  lines.  Each  of  the  fifteen  data  sections  con¬ 
tain  a  specific  number  of  variables  or  arrays  depending  upon  the  section  require¬ 
ments.  Each  variable  occupies  a  20  column  field  comprised  of  a  10  column  alph¬ 
anumeric  descriptor  space  and  a  10  column  integer  or  floating  point  data  space. 
The  descriptor  space  is  provided  for  the  user's  convenience,  and  is  usually  used 
to  enter  the  variable  name.  The  descriptor  is  typically  left- justified  with  any 
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remaining  spaces  filled  with  blanks.  The  data  space  Is  always  located  to  the 
right  of  Its  descriptor  space,  and  should  be  right-justified  in  the  designated 
f ormat * 

The  variable  fields  are  arranged  within  the  data  sections  in  records 
(or  lines)  of  from  one  to  four  fields  per  record.  The  individual  section 
requirements  for  field  arrangement  are  discussed  in  each  of  the  following  sub¬ 
sections. 


The  GESS  program  performs  a  series  of  checks  on  the  input  variables 
to  ensure  that  the  data  given  will  be  compatible  with  the  program  logic.  When 
problems  with  the  given  data  are  determined,  descriptive  error  messages  are 
printed.  These  messages  will  describe  the  problem  as  either  a  ’'Warning''  or  a 
"Fatal  Error".  Warning  messages  allow  the  program  to  continue,  but  fatal  errors 
will  halt  execution,  and  must  be  corrected  and  resubmitted.  Both  types  of  mess¬ 
ages  include  a  description  of  the  problem  encountered. 


4.3.1  Program  Controls.  The  Program  Control  section  allows  the  user  to  con¬ 
trol  the  execution  of  the  simulation.  Variable  fields  are  arranged  as  shown 
in  Figure  4-1. 


Figure  4-1.  Program  Control  Variable  Fields 


These  variables  are  described  in  Table  4-1. 
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Variable 

DOUHAT 

HEADER 

T START 

TSTOP 

ESTOP 


Table  4-1.  Program  Cbntrol  Variable  Descriptions 


Definition  (Legal  Values)  Format 

Input  delimiter  variable  A5 

("START”-beginning  of  new  input  deck) 

( "STOP" -end  of  job) 

((blank  card] -same  as  ”STOP“) 


Simulation  identifier  3(8A10) 

(any  relevant  3  lines  of  alphanumerics) 


Time  at  which  simulation  starts 
(0. 0-normal  start) 

(>0. 0-restart) 


Time  at  which  simulation  stops 
(0.0-stop  at  event  ESTOP) 

(>0. 0-stop  at  time  TSTOP) 


Event  at  which  simulation  stops 
(0-stops  a  time  TSTOP) 

(1-catapult  1  ignition) 
(2-catapult  2  ignition) 
(3-catapult  1  separation) 
(4-catapult  2  separation) 

(5-rail  separation) 

(6-rocket  1  ignition) 

(7-rocket  2  ignition) 

(8-rocket  3  Ignition) 

(9-rocket  4  ignition) 

(10-rocket  5  ignition) 

(11-rocket  6  ignition) 

(12-rocket  1  burnout) 

(13-rocket  2  burnout) 

( 14-rocket  3  burnout) 

(15-rocket  4  burnout) 

(16-rocket  5  burnout) 

(17-rocket  6  burnout) 

(18-drogue  chute  1  deployment) 
(19-drogue  chute  1  line  stretch) 

(10-drogue  chute  1  full  inflation) 

(21-drogue  chute  2  deployment) 
(22-drogue  chute  2  line  stretch) 

(23-drogue  chute  2  full  inflation) 

(24-recovery  chute  deployment) 
(25-recovery  chute  line  stretch) 
(26-recovery  chute  full  inflation) 
(27-peak  trajectory) 
(28-seat/occupant  separation) 
(29-seat/occupant  impact) 
(30-occupant  impact) 


F10.4 


F10.4 


12 


(  Cbntinued) 
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Table  4 

Variable 

ESTOP  (Gbnlxnued) 


IRESTRT 


I UNITS 

ISEATTR 

ISOSEP 

IPLQT 

IDR1FLG 


1.  Program  Control  Variable  Descriptions  (Qontinued) 


Definition  (Legal  Values) 


Format 


(31-seat  Impact) 

(32-aircraft  impact) 

(33-DART  right  line  start) 

(34-DART  left  line  start) 

(35-DART  right  line  broken) 

(36-DART  left  line  broken) 

(37-seat  separation  from  aircraft) 

(38-seat  first  motion) 

Instantaneous  restart  flag,  to  12 

create  a  data  base  of  all  trajectory 
information  at  a  specified  time  from 
which  multiple  continuations  can  be 
made 

(0-restart  file  not  created) 

(1-not  Implemented) 


Flag  controlling  units  of  measure  12 

(0-metric  units) 

( 1-English  units) 


Flag  controlling  seat  alone  trajectory  12 

(O-trajectory  not  generated) 

(1-trajectory  generated) 


Flag  controlling  seat/occupant  separation  12 
(O-separatlcn  suppressed) 

(1-separation  based  on  time) 

(2-separation  based  on  force) 


Flag  controlling  plot  file  generation  12 

(0-plot  file  not  generated) 

(1-plot  file  generated) 


Flag  controlling  Dynamic  Response  12 

Index  Calculation 
(0-DRI  not  generated) 

(1-DRI  generated) 


Units 


N/A 


N/A 


N/A 


N/A 


N/A 


N/A 
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4.3.2  Report  Flags.  The  Report  Flags  section  allows  the  user  to  select  the 
output  reports  that  will  be  required  for  the  specific  analyses  to  be  conducted. 
Variable  fields  are  arranged  as  shown  in  Figure  4-2. 


IREPTS(l), 

IREPTS(2) , 

IREPTS(3) , 

IREPTS(4) 

IREPTS(5) , 

IREPTS(6) , 

IREPTS(7) , 

• 

IREPTS(8) 

IREPTS(29) , 

IREPTS(30) , 

• 

• 

IR£PTS(31) 

Figure  4-2.  Report  Flag  Variable  Fields 
These  variables  are  described  in  Table  4-2. 

Table  4-2.  Report  Flags  Variable  Descriptions 

Definition  (Legal  Values)  Format  Units 

Flag  controlling  report  i  generator  12  N/A 

(0-report  i  not  generated) 

(1-report  i  generated) 

Refer  to  Table  5-1  for  report  titles. 

4.3.3  Integration  Time  Steps.  The  Integration  Time  Steps  section  allows  the 
user  to  adjust  the  time  increments  between  trajectory  phases.  The  precision  of 
the  trajectory  is  Inversely  affected  by  the  size  of  the  time  steps;  hence,  a 
larger  time  step,  although  more  economical  in  terms  of  computer  time,  can  result 
in  greater  inaccuracies  in  the  integration  calculations.  The  simulation  has 
been  categorized  into  three  phases: 


Variable 

IREPTS(i) 
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Phase  I  -  from  Time  ■  0  to  Seat/ Aircraft  Separation 
Phase  II  -  from  S/A  Separation  to  Seat/ Occupant  Separation 
Phase  III  -  from  S/0  Separation  to  Occupant  Impact 


It  is  recommended  that  the  earlier  phases  have  very  small  time  incre¬ 
ments,  since  this  period  is  characterized  by  very  complex  forces  and  moments 
while  later  phases  have  larger  time  increments  to  conserve  computer  time.  Thi 
variable  fields  in  this  section  are  arranged  as  shown  in  Figure  4-3. 


PI1, 

PI2, 

PI3 

DTPHAS1, 

DTPHAS2, 

DTPHAS  3 

Figure  4-3.  Integration  Time  Step  Variable  Fields 


These  variables  are  described  in  Table  4-3. 


Table  4-3.  Integration  Time  Step 

Variable 

Descriptions 

Variable 

Definition  (Legal  Values) 

Format 

Units 

PI1 

Print  frequency  -  Phase  I 

14 

N/A 

PI2 

Print  frequency  -  Phase  2 

14 

N/A 

PI3 

Print  frequency  -  Phase  3 
(Pli  >  0) 

14 

N/A 

DTPHAS l 

Integration  time  step  -  Phase 

1 

F10.4 

sec 

DTPHAS 2 

Integration  time  step  -  Phase 

2 

F10.4 

sec 

DTPHAS 3 

Integration  time  step  -  Phase 
(DTPHASi  >  0) 

3 

F10.4 

sec 
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4.3.4  Initial  Conditions 

The  Initial  Conditions  section  allows  the  user  to  create  the  precise 
escape  scenario  required  for  analysis.  Variable  fields  are  arranged  as  shown 
in  Figure  4-4. 


TEMP, 

PRESSUR, 

DENSITY 

XPOS, 

YPOS, 

ZPOS 

YAW, 

PITCH, 

ROLL 

RVEL, 

QVEL, 

PVEL 

WINDX, 

WINDY, 

WINDZ 

XACVEL, 

ZACVEL, 

CKPITHT 

NPTSAAT , 

AAT(l.l), 

AAT( 1,2), 

• 

AAT( 1 , 3) 

AAT(50, 1) , 

• 

• 

AAT(50,2) , 

AAT(50,3 

NPTSLAT 

LAT(l.l), 

LAT( 1,2), 

• 

LAT( 1,3), 

LAT(50,1) , 

• 

• 

LAT(30,2) , 

LAT(50,3 

F: 

variables  an 

Lgure  4-4.  In 

s  described  in 

ltial  Cond: 

Table  4-4 

AAT( 1,4) 


See  Section 
4. 3. 4.1 


AAT(50,4/ 


LAT(  1 ,4) 


See  Section 
4.3. 4.2 


LAT(50,4) 
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Table  4-4.  Initial  Conditions  Variable  Descriptions 


Variable 

Definition  (Legal  Values) 

Format 

Units 

Engl.(metr .) 

TEMP 

Air  temperature 

F10.4 

O  O 

F(  C) 

PRESSUR 

Barometric  pressure(> 0.0) 

F10.4 

millibars 

*> 

DENSITY 

Atmospheric  density  or  mass/unit 
volume  (>0.0) 

(0.0-effective  density  calculated  from 
TEMP) 

F10.4 

slugs/ftJ 

XPOS 

Initial  aircraft  X-position  in  EFCS 

F10.4 

ft(m) 

YPOS 

Initial  aircraft  Y-position  in  EFCS 

F10.4 

f  t(m) 

ZPOS 

Initial  aircraft  Z-position  in  EFCS 

F10.4 

ft(m) 

YAW 

Initial  aircraft  yaw  in  EFCS 

F10.4 

deg 

PITCH 

Initial  aircraft  pitch  in  EFCS 

F10.4 

deg 

ROLL 

Initial  aircraft  roll  in  EFCS 

F10.4 

deg 

RVEL 

Initial  angular  velocity  about 
aircraft  Z-Axis,  or  yaw  velocity 

F10.4 

deg/ sec 

QVEL 

Initial  angular  velocity  about 
aircraft  Y-Axis,  or  pitch  velocity 

F10.4 

deg/ sec 

PVEL 

Initial  angular  velocity  about 
aircraft  X-Axis,  or  roll  velocity 

F10.4 

deg/ sec 

WINDX 

Wind  X-velocity  in  EFCS 

F10.4 

ft/sec 
(m/ sec) 

WINDY 

Wind  Y-velocity  in  EFCS 

F10.4 

ft/ sec 
(m/ sec) 

WINDZ 

Wind  Z-velocity  in  EFCS 

F10.4 

ft/sec 

(m/sec) 

XACVEL 

Initial  aircraft  X-velocity  in  EFCS 
(>0.0) 

F10.4 

ft/sec 
(m/ sec) 

ZACVEL 

Initial  vertical  rate  of  climb 
in  EFCS  or  sink  rate 

F10.4 

ft/ sec 
(m/ sec) 

CKPITHT 

Aircraft  cockpit  screen  height  in  ACS  FI0.4 

or  protected  height  of  seat/occupant  (>0.0) 

ft/sec 
(m/ sec) 

NPTSAAT 

Number  of  points  in  aircraft  angular 
acceleration  table  (0-30) 

13 

N/A 

AAT(i.j) 

Aircraft  angular  acceleration  table  array 
(See  Section  4.3.4. 1) 

i(4F10. 

4) 

sec .deg / 
sec  ^ 

NPTSLAT 

Number  of  points  in  aircraft  linear 
acceleration  table  (0-30) 

13 

N/A 

LAT(i.j) 

Aircraft  linear  acceleration  table  array 

i(4F10. 

4) 

sec.ft/sec^ 

(sec,m/sec2) 
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4.3.4. 1.  If  NPTSAAT  ■  0,  no  AAT  values  are  given.  When  1  <  NPTSAAT  <  50, 

NPTSAAT  lines  of  AAT  values  muse  be  entered.  Each  line  of  AAT  data  consists 
of  the  time  at  which  that  line  of  data  is  initialized  (AAT(i,l)),  followed  by 
the  X,  Y,  and  Z  values  of  the  angular  acceleration  (AAT(i,2),  AAT(i ,3) ,  and 
AAT(i,4),  respectively).  The  CESS  program  linearly  interpolates  the  values 
between  points  based  on  the  time  value.  Hence,  a  typical  aircraft  angular 
acceleration  table  could  be  set  up  as  follows: 


NPTSAAT 

5 

AAT( 1,1) 

0.000 

AAT( 1 ,2) 

0.0 

AAT( 1 , 3)  5.0 

AAT( 1 , 4)  0.0 

AAT( 2,1) 

2.999 

AAT(2,2) 

0.0 

AAT( 2, 3)  5.0 

AAT( 2,4)  0.0 

AAT(3,1) 

3.000 

AAT(3,2) 

0.0 

AAT( 3 , 3)  -5.0 

AATC3.4)  0.0 

AAT(4,1) 

5.999 

AAT(4,2) 

0.0 

AAT(4,3)  -5.0 

AAT(4,4)  0.0 

AAT( 5,1) 

6.000 

AAT(5,2) 

0.0 

AAT(5, 3)  0.0 

AATC5.4)  0.0 

4.  3. 4. 2. 

The 

LAT  array 

is  set 

up  identically 

to  the  AAT  array,  based  on 

NPTSLAT. 

as  follows 

Hence , 

• 

• 

a  typical  aircraft 

linear  acceleration  table  could  be  set  up 

NPTSLAT  2 

LAT(l,l)  0.000  LAT( 1 , 2)  0.0  LAT(1,3)  0.0  LAT(1,2)  -32.2 

LAT( 2, i)  30.000  LAT(2,1)  0.0  LAT(2,3)  0.0  LAT(2,1)  -32.2 

4.3.5  Seat  Alone  Descriptors.  The  Seat  Alone  section  allows  the  user  to 
fully  describe  all  aspects  of  the  various  escape  system  ejection  seats.  Vari¬ 
able  fields  are  arranged  as  shorn  in  Figure  4-5. 
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XPOSSRP, 

YPOSSRP, 

ZPOSSRP 

XCGSA, 

YCGSA, 

ZCGSA 

IXXSA, 

IXYSA, 

IXZSA 

IYYSA, 

IYZSA, 

IZZSA 

PHISA, 

PSISA, 

THESA 

AREAS A, 

HGHTSA, 

WGHTSA 

XPOSBOT, 

YPOSBOT, 

ZPOSBOT 

XPOSSCS, 

YPOSSCS, 

ZPOSSCS 

XPOSSRP 

YPOSSRP 

ZPOSSRP 


X  position  of  seat  alone  aerodynamics  reference  point 
Y  position  of  seat  alone  aerodynamics  reference  point 
Z  position  of  seat  alone  aerodynamics  reference  point 


F10.4 
F 10. 4 
F10.4 


ft(m) 
ft(m) 
f  t(m) 


XCGSA 

YCGSA 

ZCGSA 


X  position  of  seat  alone  center  of  gravity  in  SCS 
Y  position  of  seat  alone  center  of  gravity  in  SCS 
Z  position  of  seat  alone  center  of  gravity  in  SCS 


F10.4 

F10.4 

F10.4 


f  t(m) 
f  t(m) 
ft(m) 


IXXSA 


Seat  alone  moment  of  inertia  about  X-axis  (>0.0) 


F  10.4 


IXYSA 


IXZSA 


IYYSA 


Seat  alone  moment  of  inertia  about  X-Y  translational 
axis  (>0.0) 

Seat  alone  moment  of  Inertia  about  X-Z  translational 
axis  (>0.0) 

Seat  alone  moment  of  inertia  about  Y-axis  (>0.0) 


F10.4 


F10.4 


F10.4 


slug-f t  ^ 

(kg-m2) 
slug-f t2 
(kg-m2) 
slug-f t 2 
(kg-m2) 


slug-f t 2 


IYZSA 


IZZSA 


Seat  alone  moment  of  inertia  about  Y-Z  translational 
axis  (>0.0) 

Seat  alone  moment  of  Inertia  about  Z-axis  (>0.0) 


F  10.4 


F 10. 4 


(kg-m2) 
slug-f t 2 
(kg-m2) 
slug-f t 2 
(kg-m2) 


PHISA 

PSISA 

THESA 


Angular  position  of  seat  alone  about  ACS  X-axis 
Angular  position  of  seat  alone  about  ACS  Y-axis 
Angular  position  of  seat  alone  about  ACS  Z-axls 


F  10.4 
F10.4 
F10.4 


AREASA  Aerodynamic  reference  area  of  seat  alone  (>0.0) 


F10.4  ft2(m2) 


HGHTSA  Height  or  vertical  length  of  seat  alone  (>0.0) 


F 10.4  f t(m) 


WGHTSA  Weight  of  seat  alone  after  occupant  separation  (>0.0)  F10.4  lbs(kg) 


XPOSBOT  X-position  of  seat  bottom  in  RCS 
YPOSBOT  Y-position  of  seat  bottom  in  RCS 
ZPOSBOT  Z-position  of  seat  bottom  in  RCS 


F10.4 

F10.4 

F10.4 


ft(m) 
f  t(m) 
f  t(m) 


XPOSSCS  X-position  of  SCS  origin  in  RCS 
YPOSSCS  Y-position  of  SCS  origin  in  RCS 
ZPOSSCS  Z-posltlon  of  SCS  origin  in  RCS 


F10.4 

F10.4 

F10.4 


f  t(m) 
f  t(m) 
f  t(m) 
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4.3.6  Seat/Occupant  Descriptors.  The  Seat/Occupant  section  allows  the  user 
to  fully  describe  all  aspects  of  the  seat/occupant  combination.  Variable  fields 
are  arranged  as  shown  in  Figure  4-6. 


XCGSO, 

YCGSO, 

ZCGSO 

IXXSO, 

IXYSO, 

IXZSO 

IYYSO, 

IYZSO, 

IZZSO 

AREASO, 

WGHTSO 

Figure  4-6. 

Seat/ Occupant 

Variable  Fields 

These  variables  are  described  in  Table  4-6. 


Table  4-6.  Seat/ Occupant  Variable  Descriptions 


Variable  Definition  (Legal  Values) 


Format  Units 


XCGSO 

YCGSO 

ZCGSO 

IXXSO 

IXYSO 

IXZSO 

IYYSO 

IYZS0 

1ZZSO 


X-position  of  seat/occupant  center  of 
Y-position  of  seat/occupant  center  of 
Z-positlon  of  seat/occupant  center  of 


Seat/occupant  moment  of  inertia  about 

Seat/occupant  moment  of  inertia  about 
axis  (>0.0) 

Seat/occupant  moment  of  inertia  about 
axis  (>0.0) 

Seat/occupant  moment  of  inertia  about 

Seat/occupant  moment  of  inertia  about 
axis  (>0.0) 

Seat/occupant  moment  of  inertia  about 


gravity  in  SCS 

F10.4 

f  t(m) 

gravity  in  SCS 

F10.4 

ft(m) 

gravity  in  SCS 

F10.4 

ft(m) 

X-axis  (>0.0) 

F10.4 

slug-f t2 

(kg-m2) 

X-Y  translational 

F  10.4 

slug-ft  2 
(kg-m2) 

X-Z  translational 

F10.4 

slug-ft2 

(kg-m2) 

Y-axis  (>0.0) 

F10.4 

slug-ft2 

(kg-m2) 

Y-Z  translational 

F10.4 

slug-ft2 

(kg-m2) 

Z-axis  (>0.0) 

F10.4 

slug-ft2 

(kg-m2) 

AREASO  Aerodynamic  reference  area  of  seat/occupant  (>0.0) 
w'GHTSO  Weight  of  seat/occupant  combination  (>0.0) 


F10.4  ft2(m2) 

F10.4  lbs(kg) 
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4.3.7  Occupant  Alone  Descriptors.  The  Occupant  Alone  section  allows  the  user 
to  fully  describe  the  basic  features  of  the  occupant  before  and  after  separation 
from  the  ejection  seat.  Variable  fields  are  arranged  as  shown  in  Figure  4-7. 


IXXOA, 

IXOYA, 

IXZOA 

IYYOA, 

IYZOA, 

IZZOA 

AREAQAA, 

WGHTQAB, 

WGHTQAA 

SOSEP, 

DMPGC 

Figure  4-7.  Occupant  Alone  Variable  Fields 


These  variables  are  described  In  Table  4-7. 


Table  4-7.  Occupant  Alone  Variable  Descriptions 


Variable 

Definition  (Legal  Values) 

Format 

Units 

IXXOA* 

Occupant  alone  moment  of  inertia 

about 

X-axis 

(>0.0) 

F10.4 

slug-ft2 

(kg-m2) 

IXYOA* 

Occupant  alone  moment  of  inertia 
translational  axis  (>0.0) 

about 

X-Y 

F10.4 

slug-ft2 

(kg-m2) 

IXZOA* 

Occupant  alone  moment  of  inertia 
translational  axis  (>0.0) 

about 

X-Z 

F10.4 

slug-ft2 

(kg-m2) 

IYYOA* 

Occupant  alone  moment  of  Inertia 

about 

Y-axis 

(>0.0) 

F10.4 

slug-ft2 

(kg-m2) 

IYZOA* 

Occupant  alone  moment  of  inertia 
translational  axis  (>0.0) 

about 

Y-Z 

F10.4 

slug-ft 2 
(kg-m2) 

IZZOA* 

Occupant  alone  moment  of  inertia 

about 

Z-axis 

(>0.0) 

F10.4 

slug-ft2 

(kg-m2) 

(continued) 


Input  required,  although  not  currently  used  in  program. 
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Table  4-7.  OccupaaC  Alone  Variable  Descriptions 
(Continued) 


AREAQAA 

Aerodynamic  reference  area  of  occupant  alone 
after  separation  from  seat  (>0.0) 

F10.4 

ft^(m^) 

WGHTQAB 

Weight  of  occupant  alone  before  separation 
from  seat  (>0.0)  see  Section  4. 3. 7.1 

F10.4 

lbs(kg) 

WGHTQAA 

Weight  of  occupant  alone  after  separation 
from  seat  (>0.0)  see  Section  4. 3. 7.1 

F10.4 

lbs(kg) 

SOSEP 

Seat/occupant  separation  parameter  (>0.0) 
see  Section  4. 3. 7. 2 

F10.4 

sec; lbs 
(sec ; kg) 

DMPGC 

Aerodynamic  angular  damping  coefficient 
(0.0<DHPGC< 1.0) 

F10.4 

N/A 

4.3. 7.1.  After  seat/occupant  separation,  the  weight  of  the  occupant  alone  in¬ 
cludes  the  weight  of  the  parachute  pack  and  all  emergency  gear.  Before  separa¬ 
tion,  this  extra  weight  is  considered  part  of  the  seat. 

4. 3. 7. 2.  SOSEF  determines  the  point  at  which  seat/occupant  separation  occurs. 
It  does  this  conditionally,  depending  upon  the  means  being  used  to  determine 
separation  dictated  by  ISOSEF  in  Section  1.  Hence,  if  IS0SEP»1,  then  SOSEP  will 
be  the  time  at  which  separation  occurs.  Alternately,  if  IS0SEP*2,  then  SOSEP 
will  be  the  recovery  chute  force  required  to  produce  separation. 
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4.3.8  Rail  Data.  The  Rail  Data  section  allows  the  user  to  fully  describe 
all  aspects  of  the  seat/aircraft  rail  system.  Variable  fields  are  arranged  as 
shown  In  Figure  4-8. 


RAILNTH, 

RAILANG 

ISTRL, 

NSLBKS 

KXSB, 

KYSB, 

MU  SB 

YKTOR, 

XPOSRRE , 

YPOSRRE, 

ZPOSRRE 

XPOSSLRE, 

YPOSLRE, 

ZPOSLRE 

XPOSSB(l), 

YPOSSB( 1), 

ZPOSSB( 1) 

• 

• 

• 

See  Section 

• 

• 

• 

_  4.3.8. 1 

• 

XP0SSB( 6) , 

• 

YP0SSB( 6) 

• 

ZPOSSB( 6) 
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Table  4-8.  Rail  Data  Variable  Descriptions 


Variable 

Definition  (Legal  Values) 

Format 

Units 

RAILNTH 

Length  of  rails  (>0.0) 

F10.4 

f  t(m) 

RAILANG 

Pitch  angle  of  rails  WRT  ACS  X-Axis 

F10.4 

deg 

ISTRL 

Slider  block  flag 

(O-slipper/sllder  blocks  are  on  seat) 

( 1-slipper/ slider  blocks  are  on  aircraft) 

13 

N/A 

NSLBKS 

Number  of  slipper/sllder  blocks  (0-6) 

13 

N/A 

KXSB 

Effective  rail  spring  constant 
along  RCS  X-axis(>0.0) 

F10.6 

lbs/ft 

(N/m) 

KYSB 

Effective  rail  spring  constant 
along  RCS  Y-axis ( >0.0) 

F10.6 

lbs/ft 

(N/m) 

MUSB 

Rail  to  slider  block  coefficient  of 
friction  (>0.0) 

F10.6 

N/A 

YKTOR 

Effective  torsional  spring  constant 
along  RCS  Y-axis(>0.0) 

F10.6 

ft- lbs/deg 
(N-m/deg) 

XPOSRRE 

Right  rail  attachment  point  along  ACS  X-axis 

F10.4 

ft(m) 

YPOSRRE 

Right  rail  attachment  point  along  ACS  Y-axis 

F10.4 

f  t(m) 

ZPOSRRE 

Right  rail  attachment  point  along  ACS  Z-axis 

F10.4 

ft(m) 

XPOSLRE 

Left  rail  attachment  point  along  ACS  X-axis 

F10.4 

ft(m) 

YPOSLRE 

Left  rail  attachment  point  along  ACS  Y-axis 

F10.4 

ft(m) 

ZPOSLRE 

Left  rail  attachment  point  along  ACS  Z-axis 

F10.4 

ft(m) 

XP0SSB( 1) 

Sllppe^/Sllder  block  1  X-position  in  SCS 

F10.4 

ft(m) 

YPOSSB(l) 

Slipper/ Slider  block  i  Y-positlon  in  SCS 

F10.4 

f  t(m) 

ZPOSSB(  1) 

Slipper/ Slider  block  i  Z-position  in  SCS 
(see  Section  4. 3.8.1) 

F10.4 

ft(m) 

4.3.8. 1.  If  NSLBKS  »  0,  then  the  GESS  program  assumes  that  a  continuous  "rail 
within  a  rail",  with  no  slippers  or  slider  blocks,  is  being  used;  the  KPOSSB, 
YPOSSB,  and  ZFOSSB  arrays  are  not  inputted.  Otherwise,  each  of  the  NSLBKBS 
slider  blocks  Is  located  by  a  separate  line  of  XPOSSB,  YPOSSB,  and  ZPOSSB  data. 
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4.3.9  Catapult  Data.  The  Catapult  section  allows  the  user  to  fully  describe 
all  aspects  of  the  in-aircraft  catapult  seat  launching  system.  Variable  fields 
are  arranged  as  shown  in  Figure  4-9. 


INCAT 

CATLNT(l),  CATSTK(l),  TCI(l) 

XPOSAP(l),  YPOSAP(l),  ZPOSAP(l) 

See  NPTSCT( 1 ) 

Section 

4. 3. 9. I  CATHRST( 1,1,1) ,CATHRST( 1,1,2), CATHRST( 1,2,1) ,CATHRST( 1,2,2) 

See 

•  Section 

4.3. 9.2 


|CATHRST( 1,23,1) ,CATHRST( 1,25,2)  _ | 

CATLNT(2),  CATSTK( 2) ,  TCI(2) 

XPOSAP(2) ,  YP0SAP(2) ,  ZP0SAP(2) 

NPTSCT(2) 

See 

Section  CATHRST( 2,1,1), CATHRST( 2,1,2), CATHRST( 2 ,2,1), CATHRST( 2,2,2) 

4. 3.9.1 

• 

See 

.  Section 

4. 3.9. 2 


ICATHRST( 2, 25 , 1) ,CATHRST( 2 , 25 , 2) 
ITUBEND 

KTUBE,  CTUBE,  PTUBE 

MUTUBE,  EXTLNGT 


See 

Section 
4.3. 9. 3 
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Table  4-9.  Catapult  Variable  Descriptions 


Variable 

Definition  (Legal  Values) 

Format 

Units 

INCAT 

ftimber  of  catapults  (0-2) 

13 

N/A 

CATLNT(i) 

Length  of  catapult  1  tube  (>0.0) 

FI  0.4 

f  t(m) 

CATSTK(i) 

Length  of  catapult  1  stroke  (>0.0) 

FI0.4 

f  t(m) 

TCI(i) 

Catapult  1  time  of  Ignition  (>0.0) 

F10.4 

sec 

XPOSAP(i) 

YPOSAP(i) 

ZPOSAP(i) 

Catapult  1  attachment  point  on  SCS  X-axis 
Catapult  1  attachment  point  on  SCS  Y-axls 
Catapult  1  attachment  point  on  SCS  Z-axis 

F10.4 

F10.4 

F10.4 

f  t(m) 
f  t(m) 
f  t(m) 

NPTSCT(i) 

Number  of  points  in  catapult  1  thrust  table 

13 

N/A 

CATHRST(i,j  ,1) 
CATHRST(i ,j  ,2) 

Catapult  1  thrust  table  point;  time  (>0.0) 
Catapult  1  thrust  table  point;  thrust  (>0.0) 

F10.4 

F10.4 

sec 

lbs(N) 

ITUBEND 

Flag  to  simulate  catapult  tube  bending 
(O-no  tube  bending) 

( 1-tube  bending  using  default  values) 

(2-tube  bending  using  input  values) 

13 

N/A 

KTUBE 

Thbe  spring  stiffness  constant  K  (>0.0) 

F10.4 

lb/ft 

(N/m) 

CTUBE 

Tube  spring  damping  coefficient  K  (>0.0) 

F10.4 

lb- sec/ft 
(N-sec/m) 

PTUBE 

Empirical  tube  bending  constant  (>0.0) 

F10.4 

N/A 

MUTUBE 

Coefficient  of  friction  (>0.0) 

F10.4 

N/A 

EXTINCT 

External  length  of  catapult  tube  (>0.0) 

F10.4 

f  t(m) 

4.3.9. 1.  Should  INCAT«0,  then  catapults  will  not  be  simulated,  and  no  further 
data  Is  required  In  this  section.  If  INCAT* l,  then  one  catapult  will  be  simula¬ 
ted,  requiring  only  that  catapult's  data  subsection  and  the  tube  bending  subsec¬ 
tion,  If  applicable.  If  INCAT*2,  then  applicable  all  data  subsections  will  be 
required . 
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4. 3. 9. 2.  If  NPTSCT(i)  ■  0,  Chen  no  CATRRST  data  Is  required.  Otherwise,  two 
CATHRST  array  values  (time  and  thrust)  will  be  required  for  each  NPTSCT(l) 
point,  arranged  In  lines  of  four  values  each. 

4. 3. 9. 3.  If  ITUBEND’O  or  1,  then  no  further  tube  bending  data  is  required.  If 
ITUBEND-2,  then  the  five  tube  bending  variables  must  be  entered  as  Indicated. 

4.3.10  Rocket  Data.  The  Rocket  section  allows  the  user  to  fully  describe  all 
aspects  of  the  ejection  seat  booster  rockets.  Variable  fields  are  arranged  as 
shown  in  Figure  4-10. 


4.3.10.1.  If  INRKT»0,  then  rockets  will  not  be  simulated,  and  no  further  data 
is  required  in  this  section. 

4.3.10.2.  If  INRKT  >  0,  then  INRKT  represents  the  number  of  rocket  data  sub¬ 
sections  required,  each  describing  one  of  up  to  six  rockets. 

4.3.10.3.  At  least  two  rocket  thrust  table  points  must  be  given  for  each 
rocket  subsection.  Each  point  consists  of  a  time  value  and  the  thrust  at  that 
time.  The  GESS  program  linearly  interpolates  between  points  based  on  time. 

4.3.10.4.  The  rocket  thrust  angles  are  the  direction  cosines  of  the  thrust 
vector  with  respect  to  the  seat  coordinate  system.  The  thrust  in  each  direc¬ 
tion  is  computed  by  the  program  based  on  the  total  thrust  multiplied  by  the 
cosines  of  the  rocket  angles. 
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INRKT 

KKIGN(l),  RKWGHT(l),  RKBURN(l),  RKDELY(l) 

T* 

XPOSRK(l),  YPOSRK(l),  ZPOSRK(l) 

See  Secdon 

RKALPH( 1 ) ,  RKBETA( 1 ) ,  RKGAMA( 1 ) 

_4. 3. 10. 1 

RKNPTS(l) 

See  SecCion 

4.3.10.2 

RKTHRST ( 1 , 1 , 1 ) , RKTHRST ( 1 , 1 , 2 ) , RKTHRST ( 1 , 2 , 1 ) , RKTHRST( 1 , 2 , 2 ) 

See  SecCion 

• 

• 

4.3.10.3 

• 

RKTHRST ( 1,25,1) ,RKTHRST( 1,25,2)  _ 

• 

• 

ft 

RKIGN(6) ,  RKWGHT(6),  RKBOR»(6),  RKDELY(6) 

XPOSRK(6),  YPOSRK(6) ,  ZPOSRKC6) 

RKALPH( 6) ,  RKBETA( 6) ,  RKGAMA( 6) 

See  SecCion 

4.3.10.2 

RKNPTS(6) 

RKTHRST( 6 ,1 ,1 ) ,RKTHRST( 6 , 1 , 2) , RKTHRST( 6 , 2 , 1 ) , RKTHRST( 6 , 2 , 2) 

• 

See  Section 

ft 

4.3.10.3 

ft 

RKTHRST( 6 , 25 , 1) ,RKTHRST( 6 , 25 , 2)  _ 

— 

K 
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Table  4-10.  Rocket  Variable  Descriptions 
Variable  Definition  (Legal  Values)  Format  Units 

INKRT  Number  of  booster  rockets  on  seat  (0-6)  13  N/A 


RKIGN(i)  Seat  travel  on  rail  prior  to  rocket  i  F10.4  ft(m) 

ignition  (>0.0) 

RKWGHT(l)  Weight  of  fuel  in  rocket  i  (>0.0)  F10.4  lbs(kg) 

RKBURN(i)  Burn  time  for  rocket  i  (>0.0)  F10.4  sec 

RKDELY(i)  Delay  time  for  rocket  i  ignition  F10.4  sec 

after  RKIGN(i)(>0.0) 

XPOSRK(i)  Position  of  rocket  i  in  SCS  X-axis  F10.4  ft(m) 

YPOSRK(i)  Position  of  rocket  i  in  SCS  Y-axis  F10.4  ft(m) 

ZPOSRK(i)  Position  of  rocket  i  in  SCS  Z-axis  F10.4  ft(m) 


RKALPH(i)  Thrust  angle  of  rocket  i  WRT  SCS  X-axis  F10.4  deg* 
RKBETA(i)  Thrust  angle  of  rocket  i  WRT  SCS  Y-axis  F10.4  deg* 
RKGAMA(i)  Thrust  angle  of  rocket  i  WRT  SCS  Z-axis  F10.4  deg* 


RKNPTS(l)  Number  of  points  in  rocket  i  thrust  table  (2-25)  13  N/A 

RKTHRST(i,j , 1)  Rocket  i  thrust  table  point  j  time  (>0.0)  F10.4  sec 

RlCTHRST(i,J ,2)  Rocket  i  thrust  table  point  J  thrust  (>0.0)  F10.4  lbs(N) 


*  See  Section  4.3.10.4 
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4.3.11  Drogue  Chute  Data.  The  Drogue  Chute  section  allows  the  user  to  des¬ 
cribe  the  drogue  chute  system.  Variable  fields  are  arranged  as  shown  In  Figure 
4-11. 

4.3.11.1.  If  IDR0GUE»0,  then  no  further  Information  Is  required  In  this  section. 

4.3.11.2.  If  IDR0GLS*1,  then  this  subsection  Is  not  required. 

4.3.11.3.  If  IDR0GLS-0,  then  this  subsection  Is  not  required. 

4.3.11.4.  Tvo  values  are  required  for  each  line  stretch  table  point  specified 
by  NPTSDLS.  These  points  are  arranged  In  the  order  of  ascending  velocities 
with  up  to  two  points  (or  four  values)  per  line. 

4.3.11.5.  If  IFTDR01-0,  then  no  further  Information  Is  required  in  this  sub¬ 
section. 

4.3.11.6.  Two  values  are  required  for  each  line  stretch  table  point  specified 
by  NPTDROl.  These  points  are  arranged  in  the  order  of  ascending  velocities  with 
up  to  two  points  (or  four  values)  per  line. 

4.3.11.7.  If  IDROGUE  <  3,  then  no  information  is  required  in  this  subsection. 

4.3.11.8.  If  IFTDR02-0,  then  no  further  information  is  required  in  this  subsec¬ 

tion. 


4.3.11.9.  Two  values  are  required  for  each  line  stretch  table  point  specified  by 
NPTDR02.  These  points  are  arranged  in  the  order  of  ascending  velocities  with  up 
to  two  points  (or  four  values)  per  line. 


*  _»  _•  / 
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IDROGOE 

TDDPLGY, 

DISPLOY , 

DROGLL 

DRDRAG1, 

DR0GPD1 , 

POROSD1 

AREADC, 

WGHTDC, 

CDDC 

XDROGAP, 

YDROGAP, 

ZDROGAP 

IDROGLS 

DROVELX, 

DROVELY, 

DROVELZ 

NPTSDLS 

DROGLS( 1,1), 

DROGLS( 1,2), 

DROGLSC 2 

DROGLS(25,l)  ,  DROGLS(25,2)  __ 

IFTDR01 

NPTDFT1 

DROGFT1 (1,1),  DROGFT1 ( 1,2), DROGFT1 ( 2,1) .DROGFT1 ( 2 , 2~ 


DR0GFT1(25, 1) ,DR0GFT1(25,2)  _ 

VELCON 

DRORAG2  DROGPD2 ,  POROSD2 

IFTDR02 

NPTDFT2 

DROGFT2( 1,1),  DROGFT2( 1,2), DROGFT2( 2,1), DR0GFT2C 2 , 2~ 


See  Section 
"4.3.11.1 


See  Section 
"4.3.11.2 


See  Section 
"4.3.11.3 


See  Section 
"4.3.11.4 


See  Section 
"4.3.11.5  I 


See  Section 
"4.3.11.6 


See  Section 
"4.3.11.7 


See  Section 
"4.3.11.8  I 


See  Section 
"4.3.11.9 


DR0GFT2( 25 , 1 ) ,DR0GFT2( 25 , 2) 


Figure  4-11.  Drogue  Chute  Variable  Fields 


YCYTOcnew 


**.  V  *.  *.  ^ 
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These  variables  are  described  In  Table  4-11. 


Table  4-11.  Drogue  Chuce  Variable  Descriptions 


Variable 

Definition  (Legal  Values) 

Format 

Units 

IDROGUE 

Drogue  chute  control  flag 

13 

N/A 

(0-do  not  simulate  drogue  chute) 
(1-simulate  standard  single  drogue  chute) 
(2-simulate  VELCON/duplex  drogue  system) 


TDDPLOY 

Time  of  drogue  chute  deployment  (>0.0) 

F10.4 

sec 

DISPLOY 

Seat  rail  travel  required  for 
drogue  chute  deployment  (>0.0) 

F10.4 

ft(m) 

DROGLL 

Drogue  chute  line  length  (>0.0) 

F10.4 

ft(m) 

DRDRAGi 

Drogue  chute  1  drag  coefficient  (>0.0) 

F10.4 

N/A 

DROGPDi 

Drogue  chute  i  projected  diameter  (>0.0) 

F10.4 

f  t(m) 

POROSDi 

Drogue  chute  i  effective  porosity  (>0.0) 

F10.4 

N/A 

AREA  DC 

Drogue  container/slug  reference  area  (>0.0) 

F 10. 4 

ft2(m2) 

WGHTDC 

Drogue  container/slug  weight  (>0.0) 

F  10.4 

lb(kg) 

CDDC 

Drogue  container/ slug  drag  coefficient  (>0.0) 

F10.4 

N/A 

XDROGAP 

Drogue  chute  attachment  point  along  SCS  X-axis 

F10.4 

ft(m) 

YDROGAP 

Drogue  chute  attachment  point  along  SCS  Y-axis 

F  10.4 

f  t(m) 

ZDROGAP 

Drogue  chute  attachment  point  along  SCS  Z-axis 

F10.4 

ft(m) 

IDROGLS 

Drogue  chute  time  of  line  stretch  control  flag 
(0-use  slug/container  projection  velocity)* 
(1-use  input  table  data) 

13 

N/A 

DROVE LX 

Drogue  container/slug  projection  velocity  in 

SCS  X-axis 

F10.4 

ft/sec 

(m/sec) 

DROVELY 

Drogue  container/slug  projection  velocity  in 

SCS  Y-axis 

F  10.4 

ft/sec 

(m/sec) 

DROVELZ 

Drogue  container/slug  projection  velocity  in 

SCS  Z-axis 

F 10. 4 

ft/sec 

(m/sec) 

NPTSDLS  Number  of  points  In  drogue  chute  line  stretch 
table  (2-25) 


*See  Section  4.3.11.2  and  4.3.11.3 


13  N/A 

(continued) 
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4.3.12  Recovery  Chute  Data.  The  Recovery  Chute  section  allows  the  user  to 
describe  the  recovery  parachute  systea.  Variable  fields  are  arranged  as  shown 
in  Figure  4-12. 

4.3.12.1.  If  IRECOV-O,  then  the  recovery  chute  deployment  is  to  be  suppressed, 
and  no  further  information  is  required  in  this  section. 

4.3.12.2.  If  NPTSRDT-1,  then  only  the  value  of  TRDPLOY  is  entered.  If 
1<NPTSRDT«25,  then  only  the  deployment  time  table  array,  RECOVDT,  is  entered. 

4.3.12.3.  Two  values  are  required  for  each  line  stretch  table  point  specified 
by  NPTSRLS.  These  points  are  arranged  in  the  order  of  ascending  velocities, 
with  up  to  two  points  (or  four  values)  per  line. 

4.3.12.4.  If  IFTRECVK)  or  1,  then  no  further  information  is  required  by  this 
section. 

4.3.12.5.  Two  values  are  required  for  each  full  inflation  table  point  specified 
by  IFTRECV.  These  points  are  arranged  in  the  order  of  ascending  velocities , 
with  up  to  two  points  (or  four  values)  per  line. 

4.3.13  DART  Data.  The  Directional  Automatic  Realignment  of  Trajectory  (DART) 


sections  allows  the  user  to  fully  describe  all  aspects  of  the  DART  stabilization 
system.  Variable  fields  are  arranged  as  shown  in  Figure  4-13. 
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1BEC0V 


RECDRAG 

POROSR, 

YRECAP, 

ZRECEP 

CHALT2, 

GLIMIT 

XRECAP,  YRECAP,  ZRECEP 

CHALT1,  CHALT2,  GLIMIT 

TDELAY 

NPTSRDT 

TRDPLOY 

RECOVDT( 1,1), RECOVDT( 1,2), REC0VDT( 2,1), RECOVDT( 2,2) 


RECOVDT( 25 , 1 ) , RECOVDT( 2 5 , 2 ) 

NPTSRLS 

RECOVLS( 1,1) ,RECOVLS( 1 , 2) ,R£COVLS( 2,1) ,RECOVLS( 2,2) 


RECOVLS(25,l) ,RECOVLS(25,2) 

IFTRECV 

NPTSRFT  __ 

RECOVPT( 1,1), RECOVFT( 1,2), REC0VFT( 2,1), RECOVFT( 2,2) 


See  Section 
"4.3.12.1 


See  Section 
"4.3.12.2 


See  Section 
"4.3.12.3 


See  Section 
"4.3.12.4 


See  Section 
"4.3.12.5 


RECOVFT(25, 1) ,RECOVFT( 25 , 2)  _ | 


Figure  4-12.  Recovery  Chute  Variable  Fields 
These  variables  are  described  In  Table  4-12. 
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Table  4-12.  Recovery  Chute  Variable  Descriptions 


Variable 

Definition  (Legal  Values) 

Format 

Units 

I  RE  00  V 

Recovery  chute  control  flag 
(0-do  not  simulate  recovery  chute) 

( 1-slmulate  recovery  chute  with  fixed  time  delay) 
(2-simulate  recovery  chute  with  time  delay  extended 
for  exceeding  GLIMIT) 

13 

N/A 

REGOVLL 

Recovery  chute  line  length  (>0.0) 

F10.4 

f  t(m) 

RECDRAG 

Recovery  chute  drag  coefficient  (>0.0) 

F10.4 

N/A 

REOOVPD 

Recovery  chute  projected  diameter  (>0.0) 

F10.4 

ft(m) 

POROSR 

Recovery  chute  effective  porosity  (>0.0) 

F10.4 

N/A 

XRECA.P 

Recovery  chute  attachment  point  along  SCS  X-axis 

F10.4 

f  t(m) 

YRECAP 

Recovery  chute  attachment  point  along  SCS  Y-axis 

F10.4 

f  t(m) 

Z  RE  CAP 

Recovery  chute  attachment  point  along  SCS  Z-axis 

F10.4 

f  t(m) 

CHALT1 

Lower  chute  deployment  altitude  (>0.0) 

F10.4 

f  t(m) 

CHALT2 

Upper  chute  deployment  altitude  (>CHALT1) 

F10.4 

f  t(m) 

GLIMIT 

Acceleration  limit  in  SCS  Z-axis  for  chute 
deployment  delay  (>0.0) 

F10.4 

G 

TDELAY 

Recovery  chute  deployment  delay  time  between 

CHALT1  and  CHALT2  (>0.0) 

F10.4 

sec 

NPTSRDT 

Number  of  points  in  recovery  chute  deployment 
time  table  ( 1-25) 

13 

N/A 

TRDPLOY 

Time  at  which  recovery  chute  is  deployed  (>0.0) 

F10.4 

sec 

NPTSRLS 

Number  of  points  in  recovery  chute  line 
stretch  table  (2-25) 

13 

N/A 

RE00VLS(i,  1) 

Recovery  chute  line  stretch  table  point 
i  velocity  (>0.0) 

F10.4 

ft/ sec 
(m/ sec 

RE0DVLS(i,2) 

Recovery  chute  line  stretch  table  point  i  time  (>0.0) 

F10.4 

sec 

IFTRECV 

Recovery  chute  fill  time  control  flag 
(0-use  default  table  data) 

(1-use  input  table  data) 

13 

N/A 

NPTSRFT 

Number  of  points  in  recovery  chute  fill  time  table 
(2-25) 

13 

N/A 

RE(DVFT(i,  1) 

Recovery  chute  fill  time  table  point  i  velocity 
(>0.0) 

F10.4 

ft/sec 
(m/ sec 

RECDVFT(i,2) 

Recovery  chute  fill  time  table  point  i  time  (>0.0) 

F10.4 

sec 
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LDART 
DRTFRCE, 
XDRTAP(l), 
XDRTCP(l) , 
XDRTAP(2) , 
XDRTCP(2) , 


DRTSTRT , 
YDRTAP( 1 ) , 
YDRTCP(l), 
YDRTAP(2) , 
YDRTCP(2)  , 


DRSTOP 

ZDRTAP(l) 

ZDRTCP(l) 

ZDRTAP(2) 

ZDRTCP(2) 


See  Section 
'4.3.13.1 


Figure  4-13.  DART  Variable  ^xds 


These  variables  are  described  in  Table  4-13. 


Table  4-13.  DART  Variable  Descriptions 


Variable 

Definition  (Legal  Values) 

Format 

IDART 

DART  system  flag 
(0-no  DART  system) 

( 1-DART  system) 

13 

DRTFRCE 

Force  which  the  DART  lines  exert  on  the  seat 
00.0) 

F10.4 

DRTSTRT 

Distance  between  seat  and  aircraft  at  which  the 
DART  initiates  stabilization  00.0) 

F10.4 

DRTSTOP 

Distance  between  seat  and  aircraft  at  which  the 
DART  lines  break  ODRTSTRT) 

F10.4 

XDRTAP(i) 

DART  line  i  cockpit  attachment  point 
along  SCS  X-axis 

F10.4 

YDRTAP(i) 

DART  line  i  cockpit  attachment  point 
along  SCS  Y-axis 

F10.4 

ZDRTAP(i) 

DART  line  i  cockpit  attachment  point 
along  SCS  Z-axis 

F10.4 

XDRTCP(i) 

DART  line  i  confluence  point  along  SCS  X-axis 

F10.4 

YDRTCP(l) 

DART  line  1  confluence  point  along  SCS  Y-axis 

F10.4 

ZDRTCP(i) 

DART  line  i  confluence  point  along  SCS  Z-axis 

F10.4 

4.3.13.1.  If  IDAKT“0 ,  no  further  information  is  required  in  this  section 
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4.3.14  TVC  Data.  The  Thrust  Vector  Control  (TVC)  section  allows  the  user 
to  implement  and  describe  the  response  of  a  gimbaled  booster  rocket.  Variable 
fields  are  arranged  as  shown  in  Figure  4-14. 


ITVC 

MPHI, 

MPSI, 

MTHE 

See  Section  4.3.14.1 

ROLLRL, 

PITCHRL, 

SMPLRAT 

TVCDLAY, 

RKANG, 

— 

Figure  4-14.  TVC  Variable  Fields 


These  variables  are  described  in  Table  4-14. 


Table  4-14.  TVC  Variable  Descriptions 


Variable 

Definition  (Legal  Values) 

Format 

Units 

ITVC 

Thrust  vector  control  (TVC)  flag 
(0-do  not  simulate  TVC) 

( 1-simulate  TVC) 

13 

N/A 

MPHI 

TVC  coordinate  system  rotation  about  SCS 

X-axis 

F10.4 

deg 

MPSI 

TVC  coordinate  system  rotation  about  SCS 

Y-axis 

F10.4 

deg 

MTHE 

TVC  coordinate  system  rotation  about  SCS 

Z-axis 

F10.4 

deg 

ROLLRL 

Rocket  rotation  limit  about  TVCCS  X-axis 

F10.4 

deg 

PITCHRL 

Rocket  rotation  limit  about  TVCCS  Y-axis 

F10.4 

deg 

SMPLRAT 

Gyroscope  maximum  sampling  rate  (>0.0) 

F10.4 

deg/ 

sec 

TVCDLAY 

TVC  gimbaling  time  delay  after  rocket  ignition  (>0.0) 

F10.4 

sec 

RKANG 

Neutral  rocket  pitch  angle  in  TVCCS 

F10.4 

deg 

-32 
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4.3.14.1.  If  ITVC*0,  Chen  no  further  information  is  required  in  this  section. 

4. 3. IS  Dynamic  CG  Data.  The  Dynamic  Center  of  Gravity  (CG)  section  allows 
the  user  to  implement  and  describe  the  response  of  an  occupant  modeled  as  a 
single  mass  connected  to  the  seat  by  a  three-dimensional  spring-damper  system. 
Variable  fields  are  arranged  as  shown  in  Figure  4-15. 
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T l  le  4-15.  Dynamic  CG  Variable  Descriptions 


Variable  Definition  (Legal  Values)  Format  Units 


IDYNCG 

Dynamic  CG  flag 

(0-do  not  simulate  dynamic  CG) 

( 1-simulate  dynamic  CG  using  default  values) 
(2-simulate  dynamic  CG  using  input  values) 

13 

N/A 

CX 

Damping  constant  along  SCS  X-axis  (>0.0) 

F10.4 

lbs- sec/ft 
(N-sec/m) 

XSLACK 

No-response  distance  along  SCS  X-axis  (>0.0) 

F10.4 

f  t(m) 

SXF 

Spring  constant  along  SCS  positive  X-axis  (>0.0) 

F10.4 

lbs/ft 

(N/m) 

SXN 

Spring  constant  along  SCS  negative  X-axis  (>0.0) 

F10.4 

Ibs/ft 

(N/m) 

CY 

Damping  constant  along  SCS  Y-axis  (>0.0) 

F10.4 

lbs-sec/f t 
(N-sec/m) 

SY 

Spring  constant  along  SCS  Y-axis  (>0.0) 

F10.4 

lbs/ft 

(N/m) 

CZ 

Damping  constant  along  SCS  Z-Axis  (>0.0) 

F10.4 

lbs-sec/f t 
(N-sec/m) 

Z  SLACK 

No-response  distance  along  SCS  Z-axis  (>0.0) 

F10.4 

ft(m) 

SZP 

Spring  constant  along  SCS  positive  Z-axis  (>0.0) 

F10.4 

lbs/ft 

(N/m) 

ZBOT 

Absolute  value  of  bottoming  distance  along 

SCS  Z-axis  (>0.0) 

F10.4 

f  t(m) 

SZN1 

Spring  constant  below  ZBOT  (>0.0) 

F10.4 

lbs/ft 

(N/m) 

SZN2 

Spring  constant  above  ZBOT  (>SZN1) 

FiO.4 

lbs/ft 

(N/m) 

4.3.15.1.  If  IDYNCGH3  or  1,  Chen  no  further  information  is  required  in  this 


section 
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4.4  Execution  Job  Stream 

The  GESS  program  is  executed  by  submitting  to  the  computer  a  series 
of  operating  system  commands  known  as  the  "job  stream".  The  job  stream  can  be 
entered  line-by-line  interactively.  However,  it  is  usually  more  convenient  to 
prepare  a  job  stream,  or  "submit  file",  in  advance  of  execution  and  “batch" 
submit  this  file  to  the  computer  for  execution.  The  interactive  commands  to 
batch  submit  a  job  stream  are: 

GET, filename 
S  UBMIT , f 1 lename 

Another  print  site  can  be  specified  if  the  main  computer  center  printer  is  not 
desired. 


SUBMIT, f i lename ,  EI*print  site  identifier 
If  no  printed  output  is  desired : 

SUBMIT.fi lename, N 

Figure  4-18  presents  a  recommended  job  stream  for  GESS  execution  using 
the  CDC  KRONOS  operating  system.  This  job  stream  provides  for: 

.  the  loading  of  the  prepared  input  file  (GESSI) 

.  the  loading  of  the  ACT  program-generated  aerodynamic  coefficient 

tables  (AER04) 

.  the  loading  and  execution  of  the  pre-compiled  program  binary  code 
(GESSB) 

.  the  optional  creation  or  replaceme  at  of  program-generated  plot¬ 
ting  data  files  (GESST  and  GESSZ) 

.  the  creation  of  a  single,  Inclusive,  indirect  access  output  file 
containing  a  mirror  copy  of  the  original  input  file  and  all 
requested  program-generated  reports  (GESSO) 
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.  A  “day  file"  record  of  Che  submitted  job  execution  (GESSDAY) 

.  An  error  recovery  logic  thac  saves  Che  day  file  and  any  prior  out¬ 
put,  in  Che  evencualicy  of  an  execucion  error 


The  recommended  job  scream  is  incended  for  execucion  runs  generacing 
a  nominal  number  of  oucput  reports  and/or  simulacing  a  relatively  few  seconds 
of  simulated  trajectory.  Should  a  large  number  of  reports  and/or  a  long  simula¬ 
tion  period  be  desired,  the  output  (and,  possibly,  the  plotting  files)  should  be 
created  either  as  a  direct  access  file,  or  as  several  indirect  files.  Refer  to 
Reference  12  for  details  on  direct  file  creation  on  the  KRONOS  system. 


/JOB 

GESS(CB200000,T75) 
ACCOUNT (XXXXXX, YYYYYY) 
ASSIGN, MS, OUTPUT 
GET(TAPEl-KGESI) 
GET(TAPE2-AER04) 
GET(GESSB) 

MAP(OFF) 

GESSB. 

REPLACE(TAPE42=»KGEST) 
REPLACE (TAPE4 i-KGESZ ) 
GOTO, 1. 

EXIT. 

1, REWIND,*. 

SKIPF( DAYFILE ) 

COPYSBF ( TAPE 1 , OUTPUT ) 
COPYBF ( TAPE  5 , OUTPUT ) 
COP YBF( TAPE 6, OUTPUT) 
COPYBF ( TAPE  7 , OUTPUT ) 


COPYBF ( TAPE  40 .OUTPUT ) 
REWIND, OUTPUT. 

REPLACE ( OUTPUT -KGESO ) 
GOTO, 2. 

EXIT. 

2, DAYFILE, GESSDAY. 
REPLACE, GESSDAY. 

EXIT. 

/EOR 

/EOI 


Figure  4-19.  GESS  Execution  Job  Stream  for  CDC  KR0N0S/0S 
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5.0  OUTPUT  INTERPRETATION 


5.1  Report  Descriptions 

The  CESS  program  can,  at  the  user's  option,  generate  up  to  31  differ¬ 
ent  output  reports.  Each  of  these  reports,  listed  In  Table  5.1,  describes  a 
different  aspect  of  ejection  seat  escape  simulated  performance.  The  reports 
show  specific  performance  variables  in  relation  to  the  simulated  time  of  occur- 
ance,  and  also  show  the  times  at  which  specific  sequential  events  occur.  An 
example  series  of  output  reports  is  presented  in  Appendix  B. 


Report 

Number 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 


Table  5-1.  GESS  Report  Titles 

Report  Title 
Input  Validation 

Seat/Occupant  Linear  Time  History 

Seat/ Occupant  Angular  Time  History 

Occupant  Alone  Linear  Time  History 

Seat  Alone  Linear  Time  History 

Seat  Alone  Angular  Time  History 

Seat/ Occupant  Linear  Time  History  WRT  Aircraft 

Seat/Occupant  Angular  Time  History  WRT  Aircraft 

Occupant  Alone  Linear  Time  History  WRT  Aircraft 

Seat  Alone  Linear  Time  History  WRT  Aircraft 

Seat  Alone  Angular  Time  History  WRT  Aircraft 

Catapult  Forces  and  Moments 

Rocket  1  Forces  and  Moments 

Rocket  2  Forces  and  Moments 


(continued) 
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Table  5-1.  GESS  Report  Titles  (Oont'd) 

Report 

Number 

Report  Title 

15 

Rocket  3  Forces  and  Moments 

16 

Rocket  4  Forces  and  Moments 

17 

Rocket  5  Forces  and  Moments 

18 

Rocket  6  Forces  and  Moments 

19 

DART  Forces  and  Moments 

20 

Drogue  Forces  and  Moments 

21 

Recovery  Quite  Forces  and  Moments 

22 

Thrust  Vector  Oontrol  Microprocessor 

Data 

23 

Rail  Forces  And  Moments 

24 

Seat/Occupant  Aerodynamic  Forces  and 

Moments 

25 

Occupant  Alone  Aerodynamic  Forces  and 

Moments 

26 

Seat  Alone  Aerodynamic  Forces  and  Moments 

27 

Occupant  Alone  Angular  Time  History 
(not  currently  Implemented) 

28 

Aircraft  Linear  Time  History 

29 

Aircraft  Angular  Time  History 

30 

Aerodynamic  Coefficients  Time  History 

31 

Dynamic  Center  of  Gravity  Data 

5.2  Plotting  Flies 

The  user  has  the  option  of  creating  up  to  two  sequential  files  of  simu¬ 
lated  trajectory  data  which  can  be  used  for  generating  trajectory  plots  on  a 
Calcomp-type  plotter.  When  the  IPLOT  flag  is  set  equal  to  1  (see  Section  4.3.1), 
an  EFCS-based  plotting  file  is  generated  as  TAPE42,  and  an  ACS-based  plotting 
file  is  generated  as  TAPE41  (see  Section  4.4).  The  contents  of  either  of  these 


W*  •-‘■O  V  --  •- 


NADC-83111-60 

files  Is  described  in  Figure  5-1.  The  trajectory  variables  that  are  recorded 
in  each  of  these  files  are  listed  in  Table  5-2. 

The  plotting  files  have  been  specifically  designed  to  be  compatible 
with  the  Data  Reduction  and  Analysis  System  (DRAS).  This  program  has  three  pri¬ 
mary  capabilities.  First,  data  reduction  routines  are  provided,  such  as  fre¬ 
quency  filtering,  moving  average  calculation,  interpolation,  integration,  dif¬ 
ferentiation,  and  others.  Second,  a  simplified  plot-generation  instruction  set 
are  provided.  Third,  output  data  files  may  be  provided  for  generating  input 
files  for  additional  analysis. 

An  example  of  a  DRAS-generated  plot,  comparing  a  GESS-calculated  tra¬ 
jectory  with  the  trajectory  of  an  ejection  seat  during  field  performance  test¬ 
ing,  is  shown  in  Figure  5-2.  The  use  of  the  DRAS  program  in  conjunction  with 
GESS  plotting  files  is  detailed  in  Reference  15.  Users  desiring  to  access 
GESS-generated  plotting  files  with  other  analysis  routines  should  review  sub¬ 
routines  PLOTBIN  and  PLOTWAC,  as  given  in  the  Volume  II  GESS  Programmer's 
Manual . 


Record  No. 

File  Section 

1 

Header 

l 

General  Data 

2,3 

Variable  Names 

4 

Simulation  Data 

5-n* 

Trajectory  Data 

*  one  record  per  time  point 

Figure  5-1.  Plotting  File  Structure 
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Table  5-2.  TAPE42  Plotting  File  Trajectory  Variables 

TAPE42  Trajectory  Variables  are  in  the  EFCS. 

TAPE41  Trajectory  Variables  are  in  the  ACS. 


Variable  No. 


Variable  Name 


S0/0A  x-acceleration 
S0/0A  y-acceleration 
S0/0A  z-acceleration 
S0/0A  total  acceleration 
S0/0A  x-velocity 
S0/0A  y-velocity 
S0/0A  z-velocity 
S0/0A  total  velocity 
S0/0A  x-position 
S0/0A  y-position 
S0/0A  z-position 
S0/0A  roll  rate 
S0/0A  pitch  rate 
S0/0A  yaw  rate 
S0/0A  roll  angle 
S0/0A  pitch  angle 
SO/OA  yaw  angle 
SO/SA  x-acceleration 
SO/SA  y-acceleration 
SO/SA  z-acceleration 
SO/SA  total  acceleration 
SO/SA  x-velocity 
SO/SA  y-velocity 
SO/SA  z-velocity 
SO/SA  total  velocity 
SO/SA  x-position 
SO/SA  y-position 
SO/SA  z-position 
SO/SA  roll  rate 
SO/SA  pitch  rate 
SO/SA  yaw  rate 
SO/SA  roll  angle 
SO/SA  pitch  angle 
SO/SA  yaw  angle 


GESS  PROGRAM  VALIDATION  -  RUN  NO.  4-20  MARCH  1903 


■  OO  30.00  45.00  60.00  75.00  90.00  ibjS.OO  120-00  135.00  150.00 
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6.0  PERFORMANCE  ASSESSMENT 

A  series  of  simulations  has  been  performed  using  input  variables 
extracted  from  ejection  sled  test  data  obtained  at  the  Naval  Weapons  Center 
Supersonic  Naval  Ordinance  Research  Track  (NWC/SNORT)  facility  at  China  Lake, 
CA  during  the  ESCAPAC  Replacement  Test  Program.  A  summary  of  the  input 
variables  describing  these  simulations  are  listed  in  Table  6-1  while  the  com¬ 
plete  input  listings  of  these  simulations  are  presented  in  Appendix  D.  The 
sled  tests  chosen  for  these  simulations  were  selected  based  upon  dummy  size  (3 
and  98  percentile)  and  sled  downrange  velocity  (zero,  low,  and  high  speeds). 

To  enhance  the  interpretation  of  the  simulation  and  test  data  compari¬ 
sons,  a  series  of  plots  was  generated,  comparing  each  of  the  simulation  trajec¬ 
tories  to  the  photometric  trajectory  data  taken  from  the  corresponding  sled 
test.  Four  types  of  comparison  plots  were  generated  for  each  of  four  simula¬ 
tion/test  data  sets,  and  the  four  plots  of  each  type  are  presented  in  Figures  6-1 
through  6-4.  These  figures  are  discussed  in  detail  below. 

Figure  6-1  illustrates  the  seat/occupant  (S/0)  and  occupant  alone 
(0/A)  trajectories  in  the  X-Z  plane.  It  can  be  seen  that,  despite  the  wide 
variations  in  wind,  aircraft  velocity,  and  dummy  weight,  the  simulated  tra¬ 
jectories  closely  parallel  their  corresponding  test  (actual)  trajectories. 
It  should  be  noted  that  three  of  the  four  simulations  exhibited  higher  tra¬ 
jectory  apogee  than  test  apogee.  In  simulation  2,  predicted  apogee  is  approxi¬ 
mately  50  percent  higher  than  actual.  However,  the  downrange  distance  at 
impact,  in  all  cases,  shows  excellent  agreement. 

Figure  6-2  Illustrates  the  downrange  velocity  time  histories  of  the 
S/0-0/A.  Despite  a  wide  variation  in  the  initial  ejection  velocities,  and 
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Table  6-1.  Test  Data 

Used  as 

Simulation 

Input 

Simulation  Number 

1 

2 

3 

ER  A-7  Test  System 

SRT 

DVT 

DVT 

Test  Number 

PI 

3 

2 

Test  Date* 

5/20 

4/19 

3/27 

Dummy  Percentile  (Z) 

98 

3 

98 

Dummy  Weight  (lbs) 

250.7 

174.2 

255.5 

Total  Ejected  Weight  (lbs) 

414.2 

341.2 

418.5 

S/0  CG  Location,  X  Axis  (in) 

10.7 

8.58 

10.22 

S/0  CG  Location,  Z  Axis  (in) 

12.3 

11.62 

12.43 

X  Velocity  at  Ejection  (KEAS) 

0 

442 

201 

Air  Density  (slug/ft^  x  lO-^) 

212 

211 

216 

Wind  Vector  (deg  wrt  tra j .  vector) 

220 

230 

30 

Wind  Velocity  (ft/ sec) 

6.7 

7.0 

21.9 

Air  Temperature  (deg  F) 

71.5 

75 

59 

Barometric  Pressure  (in  Hg) 

27.4 

27.4 

27.3 

Time  of  Seat  First  Motion  (sec) 

0.031 

0.020 

0.026 

Time  of  Drogue  Projection  (sec) 

0.166 

0.155 

0.158 

Time  of  Drogue  Inflation  (sec) 

1.700 

0.286 

0.394 

Time  of  Booster  Rocket  Ignition  (sec) 

0.187 

0.179 

0.188 

Time  of  Booster  Rocket  Burnout  (sec) 

0.423 

0.415 

0.436 

Time  of  Parachute  Deployment  (sec) 

0.289 

1.417 

0.333 

Time  of  Parachute  Line  Stretch  (sec) 

0.697 

1.685 

0.581 

Time  of  Seat/Occupant  Separation  (sec) 

0.968 

1.920 

0.768 

Time  of  Parachute  Full  Inflation  (sec) 

*  All  tests  conducted  during  1981. 

3.520 

2.668 

1.572 

O3inut*no«  GlinuiftllON 
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FIGURE  6-1.  Seat/Occupant  -  Occupant  Alone  Altitude  (Z-Axis)  vs  Downrange  (X-Axis)  Distance 
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the  distinctive  changes  in  the  deceleration  rates,  the  simulation  curves  show 
exceptionally  close  comparisons  with  the  actual  data. 

Figure  6-3  illustrates  the  vertical  velocity  time  histories  of  the 
S/O-O/A.  Here,  all  plots  show  an  extended  vertical  acceleration  for  the  simu¬ 
lated  curve  relative  to  the  actual  curve.  This  could  possibly  be  due  to  a 
shorter  actual  rocket  thrust  period  or  a  stronger  than  anticipated  DART  braking 
force.  The  extended  acceleration  could  easily  explain  the  higher  predicted 
trajectory  apogee  tendency  noted  in  Figure  6-1. 

Finally,  Figure  6-4  illustrates  the  total  velocity  time  history  of  the 
S/O-O/A.  It  can  be  seen  that  in  all  cases  the  agreement  between  the  simu¬ 
lated  and  actual  time  histories  is  again  excellent. 

Since  "the  purpose  of  all  computation  is  insight,  rather  than  num¬ 
bers"  (10),  it  is  essential  that  the  user  know  and  understand  the  effective 
limits  to  the  use  of  the  GESS  program.  The  foregoing  brief  analysis,  while 
not  intended  as  a  complete  validation  study,  suggests  that  the  GESS  program 
provides  good  representation  of  actual  escape  system  performance.  The  model 
has  not  been  designed  to  duplicate  the  intrinsic  stochastic  variability  that 
exists  in  real  world  ejection  situations,  and  should  not  be  used  with  that 
intent.  However,  we  conclude  that  when  correctly  used  with  reasonably  accurate 
input  data,  the  GESS  program  will  provide  highly  reliable,  detailed  descrip¬ 
tions  of  escape  system  performance. 


FIGURE  6-3.  Seat/Occupant  -  Occupant  Alone  Vertical  Velocity  (Z-Axls)  vs  Time 
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FIGURE  6-4.  Seat/Occupant  -  Occupant  Alone  Total  Velocity  vs  Time 
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APPENDIX  A 

AERODYNAMIC  COEFFICIENT  TABLES 

The  Aerodynamic  Coefficient  Tables  (ACT)  program  creates  aerodynamic 
coefficient  tables  that  are  functions  of  either  two  or  three  Independent 
variables.  The  format  in  which  the  ACT  tables  are  created  is  designed  to  be 
similar  to  the  format  in  which  the  experimental  data  are  received.  A  maximum 
of  500  entries  are  allowed  in  a  table.  For  three-way  tables,  the  coefficients 
are  functions  of  the  angle  of  attack  (ALPHA),  the  angle  of  sideslip  (BETA), 
and  mach  number  (MACH).  The  two-way  tables  may  be  functions  of  BETA  and  MACH, 
BETA  and  ALPHA,  or  ALPHA  and  MACH. 

By  specifying  the  order  in  which  the  data  must  be  read  into  memory, 
and  through  the  use  of  variable  dimensioned  arrays,  the  user  can  easily  store 
and  reference  the  tables  on  random  access  files.  As  the  coefficients  are  stored 
in  consecutive  locations  in  memory,  it  is  necessary  that  the  data  be  read  as 
follows : 

8  For  three-way  data: 

C(I,J,K)-F((((ALPHA(I),BETA(  J),MACH(K)) ,1-1 ,L) , J-l ,M)  ,K-1,N) 

*  For  two-way  data: 

C(J,K)-t(((BETA(J),MACH(K)),J-l,M) ,K-1,N) 
or  C(J,I)«F(((BETA(J) ,ALPHA(I)) ,J-1,M) ,1-1, L) 
or  C(I,K)-F(((ALPHA(I),MACH(K)) ,1-1 ,L) ,K-1 ,N) 

L,  M,  and  N  are  the  number  of  ALPHA,  BETA,  and  MACH  coordinates,  respectively. 
Proper  dimensioning  of  the  coefficient  array  facilitates  printout  of  the  tables. 


A-l 
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A.l  SEQUENCE  NUMBERS 

A  sequence  number  is  assigned  to  each  table  and  is  used  to  reference 
the  disk  address  at  which  the  table  is  stored.  During  creation  runs  the  program 
assigns  consecutive  numbers  to  the  tables  as  they  are  read.  Three-way  tables 
are  assigned  sequence  numbers  1  to  20,  and  two-way  tables  are  assigned  sequence 
numbers  21  to  50.  See  Reference  13  for  further  information  on  random  files  and 
sequence  numbers. 

When  the  tables  are  extended,  the  program  searches  the  appropriate 
info  array  (see  Section  A.  3.2)  to  determine  the  sequence  number  of  the  last 
random  file  table.  Each  table  that  is  added  is  assigned  the  next  available 
sequence  number.  During  this  search,  a  check,  is  made  to  ensure  that  the  tables 
will  not  be  extended  beyond  their  limit.  To  replace  a  table,  the  user  must 
specify  the  sequence  number  of  the  table  he  wishes  to  replace. 


A.  2  INFO  ARRAYS 

Associated  with  each  table  is  an  INFO  array,  containing  information 
about  that  table.  INF03  is  a  20  by  13  array  associated  with  three-way  tables, 
and  INF02  is  a  30  by  10  array  associated  with  two-way  tables.  These  arrays  are 
updated  and  written  to  a  random  file  after  creating,  extending,  or  replacing 
tables.  To  facilitate  referencing,  the  20  by  13  INFO  array  is  assigned  disk 
address  51,  and  the  30  by  10  INFO  array  is  assigned  disk  address  52. 

During  an  extension  run,  the  proper  INFO  array  is  interrogated  to 
determine  the  number  of  tables  already  on  random  file.  The  new  tables  are 
assigned  the  next  available  sequence  numbers,  provided  there  is  room  to  extend 
them.  In  the  replacement  mode,  the  first  word  in  the  INFO  array  for  the  given 
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sequence  number  is  checked  Co  be  sure  Chat  a  Cable  wlch  chac  sequence  number 
already  exists. 


The  concents  of  the  INFO  array,  with  J  representing  Che  sequence  num¬ 
ber  of  Che  current  cable  being  processed,  are  as  shown  below. 


For  two-way  cables : 


INF02(J, 1)  - 
INF02( J, 2)  - 
INF02( J, 3)  - 
INF02( J, 4)  - 
INF02( J, 5)  - 
INF02(J,6)  - 
INF02( J, 7)  - 
INF02( J, 8)  - 
INF02( J, 9)  - 
INF02( J, 10)- 


Cable  identifier 
number  of  BETAs 
number  of  MACHs 
BETA  delta* 

MACH  delta* 
minimum  BETA 
minimum  MACH 
maximum  BETA 
maximum  MACH 
type  of  two-way  table 


Depending  on  the  type  of  two-way  table  being  processed,  either  BETA  or  MACH, 
as  used  above,  may  be  replaced  by  ALPHA. 


For  three-way  Cables: 


INF03(J, 1) 
INF03( J, 2) 
INF03( J, 3) 
INF03( J, 4) 


table  identifier 
number  of  ALPHAS 
number  of  BETAs 
number  of  MACHs 


*  delta:  increment  between  variable  points 
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INF03(J, 5)  - 
INF03(J,6)  - 
INF03( J, 7)  - 
INF03( J, 8)  - 
INF03(J, 9)  - 
INF03( J, 10)  - 
INF03(J,  11)  - 
INF03( J, 12)  - 
INF03(J, 13)  - 


ALPHA  delta* 
BETA  delta* 
MACH  delta* 
minimum  ALPHA 
minimum  BETA 
minimum  MACH 
maximum  ALPHA 
maximum  BETA 
maximum  MACH 


A.  3  COMMON  VARIABLES 

A. 3.1  Externally  Defined 

The  following  common  variables  are  user-defined  through  the  program 

input. 


ALMAX 

BEMIN 

DELMAC 

ITYPE 

MODE 

NTYPE 

ALMIN 

DELALP 

FMT 

MAMAX 

NCOEF 

BEMAX 

DELS  El' 

IDENT 

MAMIN 

NOFTAB 

The  mnemonic  definitions  of  the  above  are  presented  in  Section  4.1.8  of  this 
report. 

A. 3.2  Internally  Defined 

Table  A-l  gives  the  definitions  for  the  common  variables  which  are 
initialized  by  the  program. 


*  delta:  increment  between  variable  points 
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Table  A-l.  ACT  VARIABLES  DEFINITIONS 


IBLNK 


IC 


IERTEST 


IJ,  IK 


INF02,  INF03 


ISEQNO 


ITABNO 


I2NO 

I3NO 

MTEST 


NAL 

NBE 

NMA 


Hollerith  word  for  blanking  out  the  table  name  in  the 
INFO  array. 

The  sequence  number  of  the  current  table.  Used  in  the 
ACT  and  INPUTT  routines  when  referencing  the  INFO  array. 

Internal  test  work,  set  by  the  ERMSG  subroutine,  which 
will  alter  the  flow  of  the  program  to  compensate  for  an 
input  error. 

■  1:  ignore  present  case,  read  next  type  1  card 

*  2:  read  excess  tables  into  a  dummy  name 

*  3:  disregard  present  table,  read  next  table 
«  4:  end  program  after  present  case 

In  a  list  of  extend  run,  they  return  the  number  of  tables 
that  have  been  written  into  the  INFO  array. 

Two  dimensional  array  which  is  equivalenced  to  the  three 
dimensional  array  INF03.  The  array  INF02  (30,10)  con¬ 
tains  the  INFO  array  for  two-way  tables.  INF03  (20,13) 
contains  the  INFO  array  for  the  tables  of  three-way  data. 
Detailed  information  concerning  the  contents  of  the  INFO 
array  may  be  found  in  Secion  A. 2. 

The  sequence  number  of  the  table  to  be  replaced. 

(see  Section  A. 7) 

Position  of  table  ISEQNO  on  the  INFO  array 
-  ISEQNO-20:  for  two-way  tables 

*  ISEQNO  :  for  three-way  tables 

The  number  of  two-way  and  three-way  tables .respectively . 
Used  as  a  counter  when  printing  the  file  dictionary. 

A  control  word  used  in  the  PTOUT  and  WRINFO  routines. 

*  0:  read  the  INFO  array  from  random  file  prior  to 

processing  the  first  table  on  a  replacement  or 
extension  run 

*  1:  do  not  read  the  INFO  array  again 

*  9:  the  PTOUT  routine  only  prints  the  file  dictionary 

The  number  of  ALPHA,  BETA,  and  MACH  numbers,  respective¬ 
ly,  that  are  contained  in  the  present  table.  Computed  by 
the  program  from  the  given  maximum,  minimum,  and  range 
values  for  each  variable. 
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A.  4  SYSTEM  ROUTINES 

Table  A-2  Uses  the  routines  called  by  the  ACT  program  which  are  avail¬ 
able  on  the  CDC  system's  library.  A  more  complete  description  of  these  routines 
can  be  found  in  Reference  13. 

Table  A-2.  System  Routines  Required  by  ACT 

OPENMS  Informs  the  operating  system  that  the  mass  storage 
file  will  be  a  random  access  file. 

READMS  Transfers  data  from  mass  storage  to  central  memory. 

SECOND  Gives  accumulated  central  processor  time. 

WRITMS  Transfers  data  from  central  memory  to  mass  storage. 

A.  5  I/O  REQUIREMENTS 

The  input/output  requirements  of  the  ACT  program  require  that  the 
following  buffer  areas  be  reserved. 


TAPE10/ INPUT 

permanent  file/ 
card  reader 

card  images 

coded 

OUTPUT 

printer 

program  output 

coded 

TAPE1 

local 

aerodynamic 

random, 

tables 

binary 

It  should  be  recognized  that  TAPE l  is  a  random  file  which  may  be  saved 
as  a  permanent  file  according  to  the  procedures  explained  in  Reference  12. 

A. 6  DECK  SETUP 

The  object  code  for  this  program  is  saved  as  a  permanent  file.  The  pro¬ 
gram  may  be  referenced  under  the  filename  ACTBIN.  Loading  of  the  program  re¬ 
quires  SOK  octal  locations  of  core  memory. 
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Job  Qjntrol  for  Creation  Run 


The  following  control  cards  are  required  to  execute  ACT,  9ave  the 
tables  that  are  created,  and  store  these  tables  on  the  random  file. 


CASE1:  (If  input  data  are  from  cards) 
JOBNAME , CB 5000 0 , T 1 0 . 


ACCOUNT ,  NUMBER,  PASSWORD. 


GET.LGO-ACTBIN. 


binary  version  of  act 


SAVE, TAPE 1-AER04. 

7/  }  multi-punch 

8/ 

9 


random  aerodynamic  co¬ 
efficient  file 


data  cards 


7/  }  multi-punch 

8/ 

9 


CASE2:  (If  input  data  are  from  a  permanent  file) 
JOBNAME, CB50000,T 10 
A  CODUNT , NUMBER , PAS SWORD 
GET.LGO-ACTBIN 
GET,  CARDS 


LGO, CARDS 
SAVE, TAPE 1-AER04 


T 
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A.6.2  Job  Oontrol  for  Extension  Run 

Using  che  files  generated  In  Che  creation  run,  Che  following  control 
cards  are  necessary  Co  add  or  replace  Cables.  This  setup  can  also  be  used  to 
provide  only  a  listing  of  the  current  tables. 

°  CASE1:  (If  input  data  is  from  cards) 

JOBNAME,  CB50000,T10. 

ACCOUNT, NUMBER, PASSWORD. 

GET.LGO-ACTBIN. 

GET ,  TAPE  l*x»^R04 . 

LGO. 

REPLACE .TAPE1-AER04. 

7/ 

8/  }  multi-punch 

9 

data  cards 
6/ 

7/  }  multi-punch 

8/ 

9 

8  CASE2:  (If  input  data  is  from  a  permanent  file) 

JOBNAME,  CB50000,T10 

A  CCD UNT ,  NUMBER ,  PASSWORD 

GET.LGO-ACTBIN 

GET, TAPE 1-AER04 

GET,  CARDS 

LGO, CARDS 

REPLACE, TAPE 1-AER04 
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If  Che  REPLACE  card  la  not  included  in  Che  above  example,  any  new 
Cables  vrlccen  Co  TAPE1  will  be  lose  when  Che  program  ends.  The  REPLACE  card 
will  noC  ocherwise  affecC  program  execuCion,  and  may  remain  in  Che  deck  even 
if  no  new  Cables  are  creaCed. 


A. 7  CARD  INPUT 

The  following  are  Che  five  differenc  card  Cypes  thac  are  recognized 
by  Che  ACT  program.  MosC  runs  will  only  require  a  subseC  of  chese,  and,  unless 
specified,  Chey  muse  be  in  Che  order  shown.  All  inCeger  data  muse  be  right- 
jusCified  in  cheir  appropriace  columns , while  all  alphanumeric  daCa  muse  be  left- 
juseif led. 


Card  ty 

pe  1 

Format 

Cblumns 

Mnemonic 

Description 

15 

1-5 

ITYPE 

Type  of  data 
■  1 :  Two-way  data 
*  2:  Three-way  data 

15 

6-10 

MODE 

Type  of  run 

•  1 :  Creation  run 
■  2:  Extension  run 


■  3 :  Replacement  run 

■  4:  Short  listing 

■  5:  Long  listing 


15 

11-15 

NOFTAB 

Number  of  tables  to  be  processed 

15 

16-20 

NTYP 

Type  of  two-way  table 

-  1:  BETA  vs  MACH 

-  2:  BETA  vs  ALPHA 

-  3:  ALPHA  vs  MACH 
(Default  is  NTYP  -1) 

On  a  list  run  (MODE  ■  +  or  5),  NOFTAB  and  NTYP  do  not  have  to  be  defined 
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Card  type 

2 

Format 

Columns 

Mnemonic 

8A10 

1-80 

EWT 

Chrd  types  3, 

4,  and  5  must 

Card  type 

2 

Format 

Oolumns 

Mnemonic 

A10 

1-10 

IDENT 

15 

11-15 

NCOEF 

15 

16-20 

ISEQNO 

Card  type 

2 

when  ITYPE 

:  -  2: 

Format 

Cblumns 

Mnemonic 

F8.0 

1-8 

ALMIN 

F8.0 

9-16 

ALMAX 

F8.0 

17-24 

DELALP 

F8.0 

25-32 

BEMIN 

F8.0 

33-40 

BEMAX 

F8.0 

41-48 

DELBET 

F8.0 

49-56 

MAMIN 

F8.0 

57-64 

MAMAX 

F8.0 

65-72 

DELMAC 

Description 

The  format  under  which  the  coefficient 
will  be  read 

be  repeated,  in  sequence,  NOFTAB  times 

Description 
Table  name 

Number  of  coefficients 

Sequence  number  of  the  table  to 
be  replaced  (when  MODE  *3) 


Description 
Minimum  ALPHA 
Maximum  ALPHA 
ALPHA  Range 
Minimum  BETA 
Maximum  BETA 
BETA  Range 
Minimum  MACH 
Maximum  MACH 
MACH  Range 
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when  ITT PE 

-  1: 

Format 

Columns 

Mnemonic 

Description 

NTYP-1 

NTYP-2 

NTYP-3 

F8.0 

25-32 

BEMIN 

Min.  BETA 

Min.  BETA 

Min.  ALPHA 

F8.0 

33-40 

BEMAX 

Max.  BETA 

Max.  BETA 

Max.  ALPHA 

F8.0 

41-48 

DELBET 

BETA  Range 

BETA  Range 

ALPHA  Range 

F8.0 

49-56 

MAMIN 

Min.  MACH 

Min.  ALPHA 

Min.  MACH 

F8.0 

57-64 

MAMAX 

Max.  MACH 

Max.  ALPHA 

Max.  MACH 

F8.0 

65-72 

DELMAC 

MACH  Range 

ALPHA  Range 

MACH  Range 

0  Card  type  5 

Format  Cblumns  Description 

Specified  by  the  user  TABLE  COEFFICIENTS 

in  Card  Type  2. 

See  Section  A. 1  for  a  further  description  of  how  the  table  coefficients 
must  be  Inputted. 

A. 8  EXECUTION  ERRORS 

Types  of  errors  that  may  be  found  Include: 

1.  Number  of  tables  greater  than  50 

2.  ffcimber  of  coefficients  greater  than  500 

3.  Invalid  use  of  replace  mode 

4.  Read  excess  tables  into  dummy  file 

5.  Two-way  tables  fully  extended 

6.  Attempt  to  extend  two-way  tables  beyond  30 

7.  Three-way  tables  fully  extended 

8.  Attempt  to  extend  three-way  tables  beyond  20 

9.  Error  in  number  of  coefficients  defined 
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The  ACT  program  provides  a  method  of  continuing  execution  of  the  pro¬ 
gram,  despite  the  discovery  of  an  input  error.  The  user  is  informed  of  the  error 
that  has  been  found,  which  will  enable  him  to  correct  it  in  a  later  run.  In  cases 
where  a  serious  error  causes  the  program  to  end, the  INFO  array  for  all  previously 
created  tables  is  written  to  the  random  file. 

Uhder  a  LIST  run  (mode  ■  4  or  5),  only  card  type  1  is  recognized.  For 
all  other  runs,  all  five  card  types  are  required.  After  each  set  of  input  is 
processed,  the  next  card  type  1  is  read.  The  program  will  terminate  upon  read¬ 
ing  an  end  of  information  card. 

The  program  is  designed  to  check  for  obvious  input  errors,  and  make 
the  appropriate  adjustments  to  compensate  for  them.  The  subroutine  ERMSG  hand¬ 
les  the  printing  of  the  error  messages ,  which  indicate  the  source  of  the  error 
and  the  action  that  will  be  taken  by  the  program.  An  attempt  is  always  made  to 
save  the  INFO  array  before  terminating  the  program  due  to  an  input  error. 

Figures  A-l  and  A-2  are  provided  as  examples  of  the  inputs  needed  to 
execute  the  ACT  program. 
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1 

1  3  3 

(14F5 

.1) 

CXS 

95 

0. 

360. 

20.  0. 

1.2 

.3 

-1.0 

-0.8  -0.4  -0.2 

0.0 

0.2  0.4 

0.8 

1.0  0.8 

0.7  0.6 

0.4 

0.0 

-0.4 

-0.6  -0.8  -0.9  -1.0 

-1.0  -0.8  -0.4  -0.2 

0.0 

0.2  0.4 

0.8 

1.0  0.8 

0.7  0.6 

0.4 

0.0 

-0.4 

-0.6  -0.8  -0.9  -1.0 

-1.0 

-0.8  -0.4  -0.2 

0.0 

0.2  0.4 

0.8 

1.0  0.8 

0.7  0.6 

0.4 

0.0 

-0.4 

-0.6  -0.8  -0.9  - 

-1.0 

-1.2 

-1.0  -0.6  -0. 2 

0.2 

0.3  0.6 

0.9 

1.0  0.9 

0.7  0.6 

0.4 

0.0 

-0.4 

-0.8  -1.0  -1.2  -1.2 

-1.3 

-1.0  -0.7  -0.2 

0.0 

0.2  0.5 

1.0 

1.2  1.0 

0.8  0.6 

0.4 

0.2 

-0.4 

-0.9  -1.1  -1.2  ■ 

-1.2 

CZS 

95 

0. 

360. 

20.  0. 

1.2 

.3 

0.3 

0.0  -0.4  -0.5  - 

-0.6 

-0.7  -0.8  -0.7 

-0.4  -0.2 

0.0  0.0 

0.3 

0.4 

0.5 

0.6  0.5  0.4 

0.3 

0.3 

0.0  -0.4  -0.5  -0.6 

-0.7  -0.8  -0.7 

-0.4  -0.2 

0.0  0.0 

0.3 

0.4 

0.5 

0.6  0.5  0.4 

0.3 

0.3 

0.0  -0.4  -0.5  -0.6 

-0.7  -0.8  -0.7 

-0.4  -0.2 

0.0  0.0 

0.3 

0.4 

0.5 

0.6  0.5  0.4 

0.3 

0.4 

0.0  -0.4  -0.6  -0.6 

-0.8  -0.7  -0.4 

0.0  0.1 

0.0  0.1 

0.2 

0.4 

0.6 

0.5  0.6  0.4 

0.2 

0.1 

0.0  -0.4  -0.7  -0.8 

-0.8  -0.9  -0.8 

-0.5  -0.2 

0.0  0.2 

0.4 

0.4 

0.6 

0.7  0.6  0.5 

0.3 

CMS 

95 

0. 

360. 

20.  0. 

1.2 

.3 

-0.10 

-0.05  0.00  0.10 

0.15 

0.15  0.15 

0.10 

f 

O 

• 

o 

i 

o 

• 

o 

1 

o 

*-< 

• 

o 

1 

00 

o 

• 

o 

1 

tri 

• 

o 

-0.02 

0.0 

0.05  0.0  -0.05-0.10 

-0.10 

-0.05  0.00  0.10 

0.15 

0.15  0.15 

0.10 

0.19-0.08-0.10-0.10- 

-0. 10- 

-0.02 

0.0 

0.05  0.0  -0.05-0.10 

-0. 10 

-0.05  0.00  0.10 

0.15 

0.15  0.15 

0.10 

0.19-0.08-0.10-0.10- 

-0. 10- 

-0.02 

0.0 

0.05  0.0  -0.05- 

-0.10 

-0.15 

-0.10  0.0  0.10 

0.15 

0.17  0.13 

0.10 

0.05  0.50 

0.10  0.15 

0.10 

0.05 

0.0 

0.90  0.0  0.90-0.15 

0.15 

-0.05  0.0  0.13 

0.15 

0.19  0.18 

0.10 

0.05-0.05-0.13-0.18- 

-0. 10- 

-0.02 

0.0 

0.0  0.02-0.10-0.15 

2 


1 


(14F5.1) 

CYS 

50 

0. 

1 

45. 

5.  0. 

0.0  -0.10 

-0.15 

-0.30 

-0.35 

-0.50 

-0.60 

-0.70 

-0.85 

-1.00 

0.0  -0.10 

-0.15 

-0.30 

-0.35 

-0.50 

-0.60 

-0.70 

-0.85 

-1.00 

0.0  -0.10 

-0.15 

-0.30 

-0.35 

-0.50 

-0.60 

-0.70 

-0.85 

-1.00 

0.0  -0.10 

-0.20 

-0.40 

-0.50 

-0.60 

-0.80 

-0.90 

-1.00 

-1.20 

0.0  -0.10 

-0.30 

-0.50 

-0.60 

-0.80 

-0.90 

-1.10 

-1.20 

-1.40 

CLS 

50 

0, 

45. 

5.  ( 

). 

0.0  +0.01 

+0.02 

+0.04 

+0.05  +0.07 

+0.09 

+0.09 

+0.09 

+0.08 

0.0  +0.01 

+0.02 

+0.04 

+0.05 

+0.07 

+0.09 

+0.09 

+0.09  +0.08 

0.0  +0.01 

+0.02 

+0.04 

+0.05 

+0.07 

+0.09 

+0.09 

+0.09  +0.08 

0.0  +0.01 

+0.02 

+0.04 

+0.05 

+0.07 

+0.09 

+0.09 

+0.09 

+0.08 

0.0  +0.01 

+0.02 

+0.04 

+0.05  +0.07 

+0.09 

+0.09  +0.09 

+0.08 

1  3 

1 

1 

(14F5.1) 

CNS 

50 

23 

0. 

» 

45. 

5.  ( 

). 

0.0  -0.01 

-0.03 

-0.05 

-0.07 

-0.08 

-0.10 

-0.12 

-0.13 

-0.14 

0.0  -0.01 

-0.03 

-0.05 

-0.07 

-0.08 

-0.10 

-0.12 

-0.13 

-0.14 

0.0  -0.01 

-0.03 

-0.05 

-0.07 

-0.08 

-0.10 

-0.12 

-0.13 

-0.14 

0.0  -0.01 

-0.03 

-0.05 

-0. 06 

-0.10 

-0.12 

-0.13 

-0.14 

-0.15 

0.0  -0.01 

-0.03 

-0.05 

-0.08 

-0.10 

-0.12 

-0.13 

-0.14 

-0.15 

Figure  A-2.  ACT  Sample  Input  -  Extension  Run 
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(PROGRAM  ASSUMES  A  CONTINUOUS  ‘RAIL  WITHIN  A  RAIL 


DATE :  1  SEP  83 

RUN  DESCRIPTION:  GESS  USERS  MANUAL  -  SAMPLE  SIMULATION  INPUT/OUTPUT 

SI  I  IS- ER  A-7  SRT  04  -  350  KEAS  '  98  ^OCCUPANT  -  HIGH  SPEED  MODE 
EAT  DESCRIPTION:  STENCEL  SEAT  -  CATAPULT .ROCKET .OART .DROGUE , CHUTE  OPERATE 
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